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In Newtonian gravitational theory, a tidal Love number relates the mass multipole moment created 
by tidal forces on a spherical body to the applied tidal field. The Love number is dimensionless, 
and it encodes information about the body's internal structure. We present a relativistic theory of 
Love numbers, which applies to compact bodies with strong internal gravities; the theory extends 
and completes a recent work by Flanagan and Hinderer, which revealed that the tidal Love number 
of a neutron star can be measured by Earth-based gravitational- wave detectors. We consider a 
spherical body deformed by an external tidal field, and provide precise and meaningful definitions 
for electric-type and magnetic-type Love numbers; and these are computed for polytropic equations 
of state. The theory applies to black holes as well, and we find that the relativistic Love numbers 
of a nonrotating black hole are all zero. 



I. INTRODUCTION AND SUMMARY 

Context of this work 

The exciting prospect of using gravitational-wave de- 
tectors to measure the tidal coupling of two neutron stars 
during the inspiral phase of their orbital evolution was re- 
cently articulated by Flanagan and Hinderer [TJ [2 . The 
idea is as follows. The orbital motion of a binary sys- 
tem of neutron stars produces the emission of gravita- 
tional waves, which remove energy and angular momen- 
tum from the system. This causes the orbits to decrease 
in radius and increase in frequency, and leads to the in- 
spiraling motion of the compact bodies. Late in the in- 
spiral the gravitational waves enter the frequency band 
of the detector, and detailed features of the orbital mo- 
tion arc revealed in the shape and phasing of the wave. 
At the large orbital separations that correspond to the 
low-frequency threshold of the instrument, the tidal in- 
teraction between the bodies is negligible, and the bod- 
ies behave as point masses. As the frequency increases, 
however, the orbital separation decreases sufficiently that 
the influence of the tidal interaction becomes important. 
The bodies acquire a tidal deformation, and this affects 
their gravitational field and orbital motion; the effect is 
revealed in the shape and phasing of the gravitational 
waves. 

Flanagan and Hinderer have provided a quantitative 
analysis of this story, and they have shown that the 
tidal coupling between neutron stars is accessible to 
measurement by the current generation of Earth-based 
gravitational- wave detectors (such as Enhanced LIGO). 
This prospect is exciting, because the details of the tidal 
interaction depend on the internal structure of each body, 
and the measurement can thus reveal important infor- 
mation regarding the compactness of each body, as well 
as its equation of state; and this information is released 
cleanly, during the inspiral phase of the orbital evolution, 



well before the messy merger of the two companions. 



Newtonian theory of tidal Love numbers 

The effect of the tidal interaction on the orbital mo- 
tion and gravitational-wave signal is measured by a quan- 
tity known as the tidal Love number of each companion 
[3J. In Newtonian gravity (see, for example, Ref. [1]), 
the tidal Love number is a constant of proportionality 
between the tidal field applied to the body and the re- 
sulting multipole moment of its mass distribution. In the 
quadrupolar case, the tidal field is characterized by the 
tidal moment £ a bif) '■— ~d a bU cy ±, in which the external 
Newtonian potential f7 ex t is sourced by the companion 
body and evaluated (after differentiation with respect to 
the spatial coordinates) at the body's center-of-mass. Be- 
cause the external potential satisfies Laplace's equation 
in the body's neighborhood, the tidal-moment tensor is 
not only symmetric but also tracefree; it is a symmetric- 
tracefree (STF) tensor. 

The quadrupole moment is Q ab := J p(x a x b — 
^6 ab r 2 ) d 3 x, where p is the mass density inside the body, 
x a is a Cartesian coordinate system whose origin is at 
the center-of-mass, and r := (<5 ah x a x f> ) 1 / 2 is the dis- 
tance to the center-of-mass; the quadrupole moment is 
another STF tensor. In the absence of a tidal field the 
body would be spherical, and its quadrupole moment 
would vanish. In the presence of a (weak) tidal field, 
the quadrupole moment is proportional to the tidal field, 
and dimensional analysis requires an expression of the 
form Q a b — — |fc2-R 5 £ a 6- (We use relativistic units and 
set G = c = 1.) Here R is the body's radius, and the 
factor of | is conventional; the dimensionless constant ki 
is the tidal Love number for a quadrupolar deformation. 
Using these expressions, the Newtonian potential outside 
the body can be written as a sum of body and external 
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potentials, and we have 
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-- l -[l + 2k 2 (R/rf]£ ab (t)x a x b . (1.1) 
r 2 L J 



The first term is evidently the monopole piece of the 
potential, which depends on the body's mass M. Within 
the square brackets, the first term represents the applied 
tidal field, and the second term is the body's response, 
measured in terms of the Love number k 2 . 



In Eq. (1.1 1 the total potential was truncated to the 



leading, quadrupole order in a Taylor expansion of the 
external potential; additional terms would involve tidal 
moments of higher multipole orders, and higher powers 
of the coordinates x a . When the tidal field is a pure 
multipole of order I, Eq. (1.1 1 generalizes to 
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l + 2k l (R/ry t+1 ]£ L (t)x\ (1.2) 



Here ki is the Love number for this multipolar config- 
uration, and L := a\a 2 • • • ai is a multi-index that con- 
tains a number I of individual indices. The tidal mo- 
ment is now defined by £c(t) := — Sti/ext/G — 2)!, and 
it is symmetric and tracefree in all pairs of indices. We 
also introduced x L := x ai x°" 2 ■ ■ ■ x ai . In this general- 
ized case the l-pole moment of the mass distribution 
is the STF tensor Q L := J px^ d 3 x, where the angu- 
lar brackets indicate that all traces must be removed 
from the tensor x L ; it is related to the tidal moment 
by Ql = -[2(1 - 2)1/(2/ - l)!!]*,^ 1 ^. 



Purpose of this work 

Our purpose in this paper is to introduce a precise 
notion of tidal Love numbers in general relativity, some- 
thing that was not pursued in the original work by Flana- 
gan and Hinderer [TJH]. In fact, we provide precise defi- 
nitions for two types of tidal Love numbers: an electric- 
type Love number k e \ that has a direct analogy with 
the Newtonian Love number introduced previously, and 
a magnetic-type Love number fc mag that has no analogue 
in Newtonian gravity. Magnetic-type Love numbers were 
introduced in post-Newtonian theory in the works of 
Damour, Soffel, and Xu [5] and Favata [6]. Our defi- 
nitions apply to gravitational fields that are arbitrarily 
strong, and to (weak) tidal deformations of any multipo- 
lar order. 

Our relativistic Love numbers are defined within 
the context of linear perturbation theory, in which an 
initially-spherical body is perturbed slightly by an ap- 
plied tidal field. Our definitions are restricted to slowly- 
changing tidal fields; this means that while a tidal mo- 
ment such as £z,(t) does depend on time, to reflect the 
changes in the external distribution of matter, the de- 
pendence is sufficiently slow that the body's response 
presents only a parametric dependence upon time. This 



allows us to ignore time-derivative terms in the field equa- 
tions, because they are much smaller than the spatial- 
derivative terms. For all practical purposes the pertur- 
bation is stationary, and t appears as an adiabatic pa- 
rameter. 

Gravitational perturbations of spherically-symmetric 
bodies are described by a metric perturbation p a p that 
can be decomposed into tensorial spherical harmonics; 
each multipole can be considered separately. The com- 
plete spacetime metric is g a p — g a g + p a p, with g° a p de- 
noting the (spherically-symmetric) metric of the unper- 
turbed body. We work in the body's immediate neigh- 
borhood, and the external bodies that create the (multi- 
polar) tidal field are assumed to live outside this neigh- 
borhood. To define the relativistic Love numbers it is 
sufficient to consider the vacuum region external to the 
body, and to construct g a p in this region only; this met- 
ric will be a solution to the vacuum field equations, and 



will represent the relativistic generalization of Eq. (1.2 1 



To compute the Love numbers it is necessary to construct 
g a p in the body's interior also, and this requires the for- 
mulation of a stellar model. The external problem there- 
fore applies to any type of body, while the internal prob- 
lem refers to a specific choice of equation of state. 



External problem 

We review the external problem first. We erect a co- 
ordinate system (v,r,0,(j>) that is intimately tied to the 
behavior of light cones: The advanced-time coordinate v 
is constant on past light cones that converge toward the 
center at r = 0, r is both an areal radius and an affinc- 
parameter distance along the null generators of each light 
cone, and the angular coordinates 6 A = (9, <j>) are con- 
stant on each generator. This choice of coordinates is 
inherited from previous work on the tidal deformation of 
black holes [7]. 

In these coordinates the external metric of the unper- 
turbed body is given by ds 2 ) = —fdv 2 + 2dvdr + r 2 d£l 2 , 
in which / := 1 - 2M/r and dtt 2 := d6 2 + sin 2 0# 2 ; 
this is the Schwarzschild metric presented in Eddington- 
Finkelstein coordinates. To construct the perturbation 
we impose the light-cone gauge conditions p vr = p rr = 
Pre = Pr<t> = to ensure that the coordinates keep 
their geometrical meaning in the perturbed spacetime 
[5]. (This property makes the light-cone gauge superior 
to the popular Regge- Wheeler gauge, which does not pro- 
vide the coordinates with any geometrical meaning.) A 
perturbation of multipole order I can be decomposed into 
even-parity and odd-parity sectors, and each sector must 
be a solution to the Einstein field equations linearized 
about the Schwarzschild metric. 

The even-parity sector is generated by the electric- 
type tidal moment £l(v), an STF tensor defined in a 
quasi-Cartesian system x a related in the usual way to 
the spherical coordinates (r, 9 A ). The (21 + 1) indepen- 
dent components of this tensor can be encoded in the 
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functions £m{v), in which the azimuthal index m is an 
integer within the interval — Z < m < Z; the encoding is 
described by £ L x L = r l £ m £^Y lm (9 A ), in which Y lm 
are the usual spherical-harmonic functions. We define 
the tidal potentials 

E^{v,6 A )=Y J ^{v)Y lm {e A ), (1.3a) 

m 

^'(»/) - \ E f m(«)>l ,n (0, (L3b) 

m 

4*m a ) = jorE^Wiifl^), (i.3c) 

m 

in which and Y^g are vector and tensor spherical 
harmonics of even parity; these are defined in Sec. II. 

The odd-parity sector is generated by the magnetic- 
type tidal moment Bl[v), another STF tensor whose in- 
dependent components can be encoded (as previously) in 
the functions Bid (v) . The odd-parity tidal potentials are 



with 



B<i\v,6 A ) = ^B&\v)X l r(e A ), 



(1.4a) 



B%(v,0 A ) = X( ^Y, B ^ v ) x1 ab{0 A ), (L4b) 



in which X 1 ™ and X l ™ B are vector and tensor spherical 
harmonics of odd parity; these also are defined in Sec. II. 
There is no scalar potential B^ in the odd-parity sector. 

The metric outside any spherical body deformed by a 
tidal environment characterized by the tidal moments £l 
and Bl is calculated in Sec. III. It is given by 



9vv / 
9vr 1 j 
9vA 



(I -1)1 



r l ei (r)£<- l \ 



_ r W e4 (r)4° 
-r'+ 1 6 4 (r)^ ) , 



(i-m 

2 



3(1-1) 



9ab = r il AB 



1(1+1) 



r l+2 e 7 (r)£V 



+ l^%{r)B%. 



(1.5a) 
(1.5b) 

(1.5c) 
(1.5d) 



The radial functions are 
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-—F(-l,-l-l;-2l;2M/r), (1.7e) 
B 7 := y^—Fil + 1, Z + 1;2Z + 2; 2M/r) 
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Z + 2 



F(Z,Z + l;2Z + 2;2A-f/r). (1.7f) 



Here i? is the body's radius, and F(a, b; c; z) is the hy- 
pergeometric function. The functions A n are finite poly- 
nomials in 2M/r, while the functions B n have non- 
terminating expansions in powers of 2M/r; for selected 
values of I they can be expressed in terms of elementary 
functions such as ln(l — 2M/r) and finite polynomials 
(see Table I in Sec. III). Each one of these functions goes 
to one as r goes to infinity. And while A n is finite at 
r = 2M, we observe that B n diverges logarithmically 
when r — » 2M. 



The metric of Eqs. (1.51 is valid in a neighborhood 



of the deformed body, and it provides a definition for 
the electric-type Love numbers fc e i and the magnetic-type 
Love numbers fc mag ; these refer to the multipole order Z, 
but we suppress the use of this label to keep the notation 
clean. While the definitions seem to rely on a specific 
choice of gauge for the metric perturbation, we prove in 
Sec. Ill that our Love numbers are gauge-invariant. 

When the tidal moments are switched off the met- 
ric reduces to the Schwarzschild metric expressed in the 
light-cone coordinates (v,r, 9 A ). When the mass pa- 
rameter M is set equal to zero the metric describes 
the neighborhood of a geodesic world line in a Ricci- 
flat spacctimc. In this limit the tidal moments can 
be related to the derivatives of the Weyl tensor eval- 
uated at r = 0. According to Eqn. (1.3) of Ref. |9J, 
we have that £ L = [{I - 2)!]" 1 (Ct aita2;a3 ... a ,) STF and 



Br 



(l + l)(l-2)l]-\e aibc C 



be 



. ai ) STF , where e ab 



is the permutation symbol and the tensor components are 
listed in the quasi-Lorentzian coordinates (t := v — r, x a ); 
the STF superscript indicates that the a n indices are 
symmetrized and all traces are removed. In the space- 
time of Eq. (1.5 1 the tidal moments £l and Bl retain a 



similar relationship with the Weyl tensor, with the un- 
derstanding that the relations are now approximate and 
refer to the asymptotic behavior of the Weyl tensor for 
r > M. 



The perturbed metric of Eq. (1.51 can be compared 
with the Newtonian potential of Eq. (1.2). We define 
an effective Newtonian potential U c s by g vv =: —(1 — 
2U e s), and our expression for g vv implies that in general 
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relativity, 
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eff 



(1-1)1 



A 1 + 2k el (R/r) 2l+1 B 1 £ L {v)x 



(1.8) 

In the nonrelativistic limit, A\ and B\ are bot h ap prox- 
imately equal to unity, and we recover Eq. (1.2 1; the 
electric-type Love number k a \ reduces to the Newtonian 
number fc;. In the strong-field regime we still recognize 
the A\ term as coming from the applied tidal field, while 
the B\ term is clearly associated with the body's re- 
sponse. There is no confusion between these terms, be- 
cause the structure of A\ is that of the finite polynomial 
1 + • • • + X(2M/r) 1 , which does not contain a term of 
order (2M/r) 2l+1 ; A is a numerical factor that can be 
determined by expanding the hypergeometric function. 
Because r is geometrically well defined, we can always 
distinguish the tidal terms from the body terms in the 
metric. 

The light-cone coordinates (v, r, 9 A ) are well-behaved 
across an eventual event horizon of the perturbed space- 
time, and our formalism is capable of handling black 
holes as well as material bodies. In general, however, 



the metric of Eqs. (1.5) is not regular at the event hori- 



zon, because of the presence of the B n functions, which 
diverge logarithmically in the limit r — > 2M . To repre- 
sent a perturbed black hole the metric must be devoid of 
these terms, and this can be accomplished by assigning 
k c \ = fcmag = to a black hole. This is one of the major 
conclusions of this work: The relativistic Love numbers 
of a nonrotating black hole are all zero. This result is 
contained implicitly in Ref. [7] , but the formalism of this 
paper permits a much clearer articulation of this prop- 
erty. 



Internal problem 

To compute the relativistic Love numbers for a selected 
stellar model requires the construction of the internal 
metric (also expressed as a sum of unperturbed solution 
and linear perturbation) and its matching with the ex- 
ternal metric at the perturbed boundary of the matter 
distribution. We carry out this exercise in Sees. IV and 
V, adapting the formalism of Thorne and Campolattaro 
[TU] to our light-cone coordinates. We take the body to 
consist of a perfect fluid with a polytropic equation of 
state 



p = Kp l+l ' r 



(1.9) 



Here p is the fluid's pressure, p its proper energy density, 
K is a constant, and n is the polytropic index (another 
constant). 

Our results are presented in Figs. 1-8, and tables of 
values are provided in the Appendix. In each figure we 
plot the Love number for a selected multipolc order (from 
I = 2 to I = 5) , and for selected values of the polytropic 
index n (from n = 0.5 to n = 2.0), as a function of the 




FIG. 1: Electric-type Love numbers for I = 2, plotted as func- 
tions of the compactness parameter 2M/R. The uppermost 
curve corresponds to n = 0.5 and the stiffest equation of state. 
The lowermost curve corresponds to n = 2.0 and the softest 
equation of state. The curves in between are ordered by the 
value of n. The arrangement is the same in all other figures. 



stellar compactness parameter C :— 2M/R; this ranges 
from C = — a weak-field, Newtonian configuration — 
to C = C max , with C ma x representing the compactness of 
the maximum-mass configuration for the selected equa- 
tion of state. 

For the electric-type Love numbers we observe the fol- 
lowing features, (i) At C = we recover the Newtonian 
values for polytropes, as tabulated by Brooker and Olle 
[IT] , (ii) For a constant C, k cl decreases as the poly- 
tropic index increases; this reflects the fact that as n 
increases, the matter distribution becomes increasingly 
concentrated near the center, which inhibits the devel- 
opment of large multipole moments, (iii) For a constant 
n, fed decreases as the compactness parameter increases; 
this reflects the fact that as C increases, the strength 
of the internal gravity increases, which produces an in- 
creased resistance to tidal deformations. 

For the magnetic-type Love numbers we observe the 
following features, (i) At C — the Love numbers are 
all zero; this reflects the fact that the magnetic-type tidal 
coupling is a purely relativistic effect that has a vanishing 
Newtonian limit, (ii) For a constant C, fc mag decreases 
as the polytropic index increases; this is explained as in 
the preceding paragraph, (iii) For a constant n, fc mag 
first increases as C increases, but then it decreases af- 
ter reaching a maximum; this reflects the fact that the 
magnetic-type tidal coupling is the result of an internal 
competition: a strong field is required to produce an ef- 
fect in the first place, but it eventually causes a large 
resistance to tidal deformation. 



Damour and Nagar 

After this work was completed we witnessed the ap- 
pearance of an article by Damour and Nagar |12j in which 
almost identical work is presented. Their paper, like ours, 
is concerned with the tidal deformation of compact bod- 
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FIG. 2: Magnetic- type Love numbers for 1 — 2, plotted as 
functions of the compactness parameter 2M/R. 
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FIG. 3: Electric-type Love numbers for 1 = 3. 



FIG. 6: Magnetic-type Love numbers for 1 = 4. 
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FIG. 7: Electric-type Love numbers for I — 5. 




FIG. 4: Magnetic- type Love numbers for I — 3. 



ics in full general relativity, and presents precise defini- 
tions for electric-type and magnetic-type Love numbers. 
And their paper, like ours, presents computations of Love 
numbers for selected matter models. Their coverage of 
the parameter space is wider: Damour and Nagar ex- 
amine two types of polytropic equations of state, and 
two tabulated equations of state for realistic nuclear mat- 
ter. In addition, Damour and Nagar define and compute 
"shape Love numbers," something that we did not pursue 
in this work. 

There are superficial differences between our treat- 
ments. One concerns the choice of coordinates: Damour 




FIG. 5: Electric- type Love numbers for 1 — 4. 



FIG. 8: Magnetic-type Love numbers for I — 5. 
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and Nagar work in Schwarzschild coordinates and adopt 
the Regge- Wheeler gauge for the metric perturbation; we 
work in Eddington-Finkclstcin coordinates and the light- 
cone gauge. Another concerns notation: we adopt differ- 
ent normalization conditions for the Love numbers and 
the tidal moments. These differences are not important. 

A more significant difference concerns the conclusion 
that the tidal Love numbers of a black hole must be zero. 
In this paper we boldly proclaim this conclusion, which 
we firmly believe to be a correct interpretation of our 
results. Damour and Nagar, however, shy away from 
the conclusion, although they agree with us on the basic 
results. We do not understand the reasons behind this 
reluctance. Damour and Nagar comment on the need to 
understand "diverging diagrams that enter the compu- 
tation of interacting black holes at the 5-loop (or 5PN) 
level" before reaching a conclusion. But since the results 
presented here do not rely at all on a post-Newtonian ex- 
pansion of the field equations, the fate of 5PN terms in a 
post-Newtonian representation of interacting black holes 
seems to us to be irrelevant. We point out, also, that the 
Damour-Nagar work does not provide a very clean foun- 
dation for the tidal deformation of black holes, because 
their coordinate system is ill-behaved on the event hori- 
zon. Our light-cone coordinates were selected precisely 
because they permit a unified treatment of material bod- 
ies and black holes. 

Aside from this issue of interpretation, and as far as 
we can judge, the results presented here are in complete 
agreement with the Damour-Nagar results. The Damour- 
Nagar work was carried out in complete independence 
from us, and our work was carried out in complete inde- 
pendence from them. The near-simultaneous completion 
of our works provides evidence that the problem is in- 
teresting and timely, and the agreement is a reassuring 
confirmation that each team performed their calculations 
without error. 



Fang and Lovelace 

The deformation of a black hole produced by an ap- 
plied tidal field was previously examined by Fang and 
Lovelace [13] . who concluded that Q a b — when the 
perturbation is expressed in Regge- Wheeler gauge. Fang 
and Lovelace therefore anticipated our result that the 
quadrupolc, electric- type Love number of a black hole is 
zero. These authors, however, qualified their conclusion 
by raising doubts about the gauge invariance of the re- 
sult, and claiming that the induced quadrupole moment 
of a tidally deformed black hole is inherently ambiguous. 
We do not share these reservations. 

We first discuss the issue of gauge invariance. The 
argument advanced by Fang and Lovelace in favor of 
a gauge dependence of the tidal Love number goes as 
follows. In Newtonian theory, the coordinate transfor- 
mation r — f[l + 2x(R/f) 5 ] 1 / 2 , where \ is an arbi- 
trary constant, turns a pure tidal potential £ a bX a x b into 



[1 + 2x(R/r) 5 ]£ a f,x a x b , which appears to describe a sum 
of tidal and body potentials; the transformation shifts the 
Love number by x- Fang and Lovelace correctly dismiss 
this coordinate dependence as irrelevant in Newtonian 
theory, because r has a well-defined meaning, but they 
point out that in a relativistic context, the coordinate 
transformation could be viewed as a change of gauge. 
The implication, then, is that the relativistic Love num- 
ber can be altered by a gauge transformation. Notice 
that the argument applies to all types of compact bod- 
ies: material bodies and black holes. 

We do not accept the validity of this argument. The co- 
ordinate transformation considered by Fang and Lovelace 
is not of a type that can be associated with a gauge trans- 
formation of the perturbation theory. A gauge trans- 
formation necessarily involves coordinate displacements 
that are of the same order of magnitude as the perturba- 
tion field. But the transformation from r to r does not 
involve the perturbation at all, and represents a large 
change of the background coordinates. The new coordi- 
nate f does not share the geometrical properties of the 
original r, and one would easily be able to distinguish 
the two coordinate systems. The argument, therefore, 
does not make a case for the gauge dependence of the 
Love numbers. And in fact, the gauge invariance of k e \ 
and fc mag for all types of compact bodies (material bodies 
and black holes) is established in Sec. III. 

We next discuss the issue of ambiguity. Unlike Fang 
and Lovelace, we believe that the relativistic Love num- 
bers of compact bodies, as defined in this paper, are well 
defined and completely devoid of ambiguity. The reason 
is that the metric of Eqs. ( 1.5 1, which is presented in co- 



ordinates that have clear geometrical properties, defines 
a perfectly well-defined spacetime geometry. Given this 
spacetime, one could in principle monitor the motion of 
test masses and light rays and thereby measure its de- 
tailed features, including the mass M, the tidal moments 
£l and Bl, and the Love numbers. These measurements 
would contain no ambiguities. 

The ambiguity identified by Fang and Lovelace con- 
cerns the coupling of Q a b, the induced quadrupole mo- 
ment, to £ a bci the octupole moment of the applied tidal 
field. According to Newtonian ideas, this coupling should 
lead to a force F a = —^£ a bc Q bc acting on the compact 
body. (Once more the argument applies to all types of 
compact bodies.) Fang and Lovelace associate F a with 
P a (r), the rate of change of three-momentum contained 
within a world tube of radius r that surrounds the com- 
pact body; this is calculated by integrating the flux of 
Landau-Lifshitz energy-momentum pseudotensor across 
the world tube. They observe that the result is indeed 
proportional to £ a bc Q bc , but that the coefficient in front 
depends on r. They interpret this as a statement that the 
force is ambiguous, assign the ambiguity to Q a b, and con- 
clude that the induced quadrupole moment of a tidally- 
deformed compact body is inherently ambiguous. 

We believe that the ambiguity in P a {r) is genuine — 
the result does depend on the world tube's radius. It is 
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hasty, however, to conclude from this that F a itself is 
ambiguous, because force calculations that rely on tech- 
niques of matched asymptotic expansions |14l 115] must 
involve a limiting procedure in which both M and r are 
taken to approach zero. Although ambiguities remain in 
this procedure, they are much smaller than those claimed 
by Fang and Lovelace. At the accuracy level of our cal- 
culations, the induced quadrupole moment of a tidally- 
deformed compact body is not ambiguous. 



Suen 



An earlier determination of the induced quadrupole 
moment of a tidally deformed black hole was made by 
Suen |16j . who examined the specific case of a black 
hole perturbed by an axisymmetric ring of matter. Suen 
found that the black-hole quadrupole moment is Q ab = 
+ 4 ^M 5 £ ab , so that it gives rise to a negative Love num- 



ber, fc e 



This result contradicts our own results. 



Suen's result is wrong. The starting points of Suen's 
analysis is the perturbed metric presented in Eq. (2.6) of 
his paper. It is easy to show that while the metric does 
indeed satisfy the Einstein field equations (up to terms 
that are quadratic in the small parameter A) , it fails to be 
regular at the event horizon. The metric does not, there- 
fore, represent a perturbed black hole, and the nonzero 
result for k c \ is a consequence of this fact. The regu- 
larity of the metric perturbation p a p at r — 2M can be 
judged by examining its components in the light-cone co- 
ordinates (v, r, 9, 4>), which are regular on the event hori- 
zon. A simple calculation reveals that in Suen's notation, 
p rr = -2(2U-V)/f, where / = \-2Mjr. This is singu- 
lar at r = 2M unless 2U—V vanishes there, but Eqs. (2.7) 
of Suen's paper show instead that 2U — V — ► AM 2 in the 
limit. The perturbation is singular. 



II. TIDAL MOMENTS AND POTENTIALS 

A spherical stellar model is perturbed by an external 
tidal field characterized by the electric-type tidal mo- 
ments £l(v) and the magnetic-type tidal moments Bl(v)- 
These are symmetric-tracefree (STF) tensors, and L is a 
multi-index that contains a number I of individual in- 
dices. The tidal moments depend on v (and not on the 
spatial coordinates), but this time dependence is taken 
to be so slow that all ^-derivatives will be ignored in the 
Einstein field equations. 

We begin our discussion of tidal potentials by adopt- 
ing quasi-Cartesian coordinates x a related in the usual 
way to our spherical coordinates (r, 6 ). We write 
the transformation as x a = rfl a (9 A ), with fl a — 
[sin(9cos0,sin#sin</>, cos#] denoting the unit radial vec- 
tor. We introduce 

lab ■= S ab - Cl a Cl b (2.1) 

as the projector to the transverse space orthogonal to 
tt a , and we let fl a A := dtt a /d9 A . We note the helpful 
identities 

fl a fl A = 0, (2.2a) 

— lab^A^B — ^ab^A^Si (2.2b) 

n AB n a A n b B = 1 ah . (2.2c) 

Here £Iab — diag[l, sin 2 9] is the metric on the unit two- 
sphere, and Q AB is its inverse. We introduce Da as 
the covariant-derivative operator compatible with Qab, 
and cab as the Levi-Civita tensor on the unit two-sphere 
(with nonvanishing components eg^ = — e^g = sin#). In 
addition to Eqs. (2.2 1 we also have 



£ab 
e A a B 

D A D B n a 



tabc^A^ B^ C i 

D B D A n a = ~n a n AB . 



(2.3a) 
(2.3b) 
(2.3c) 



Organization of the paper 

In the remaining sections of this paper we present the 
details of our analysis, and describe how the results re- 
viewed previously were obtained. We begin in Sec. II 
with a discussion of tidal moments and tidal potentials, 



and motivate the definitions presented in Eqs. (1.3 1 and 



(1.4). In Sec. Ill we solve the external problem, and show 



Here and below, upper-case latin indices are raised and 
lowered with £l AB and O^, respectively. Finally, we 
note that Dc^ab — Dc^ab = 0. 

For an electric- type tidal moment £l of degree I > 2, 
the Cartesian version of the tidal potentials are defined 

by 



f(l) 
oil) 



= 2 la c lb d £ cdL _ 2 n L - 2 + lab £ {l) . 



(2.4a) 
(2.4b) 
(2.4c) 



that the metric of Eqs. (1.5 1 is a solution to the vacuum TT c n\ ■ , , . , c (i) ■ 

. LJ . Here £ K ' is a scalar potential, £ a ' is a transverse vector 



field equations linearized about the Schwarzschild metric 
In Sec. IV we formulate the internal problem for general 
stellar models, and we specialize this to polytropes in 
Sec. V. In Sec. VI we review the numerical techniques 
that were employed to generate the figures and the ta- 
bles displayed in the Appendix. 



potential, and sjQ is a transverse-tracefree tensor poten- 
tial. The angular version of the tidal potentials are 



£ {l) =£ L £l L , 

£ { A l) ■.= £Vn A = n A £ a L-itt 



L-l 



(2.5a) 
(2.5b) 



8 



£ 



(i) 



U; — ^ab^A^B — ^A^B£abL-2^ L 2 + £Iab£^ ■ 

(2.5c) 



For a magnetic-type tidal moment Bl of degree Z > 2, 
the Cartesian version of the tidal potentials are defined 

by 



B« := 



L-l 



B 



(0 

f/6 



f n p B 9 



L-2 



(2.6a) 
(2.6b) 



Here b£° is a transverse vector potential, and B a l is a 
transverse-tracefree tensor potential; there is no scalar 
potential in the magnetic case. The angular version of 
the tidal potentials are 



B 



(0 



B$n\ = n A e apq WBY 



,0 



L-l 



B 



(0 

AB 



n(0o° o b 

D ab il A il B 
b 



^B^bpq^VB q a j 



(2.7a) 



(2.7b) 



The tidal potentials can all be expressed in terms of 
(scalar, vector, and tensor) spherical harmonics. Let 
Y lm be the standard (scalar) spherical-harmonic func- 
tions. The vector and tensor harmonics of even parity 



are Y\ m := D A Y im , n AB Y l 



and Fi'g 



:= [D A D B 



\l(l + l)Q AB ]Y lm ; notice that n AB Y A ^ = by virtue 
of the eigenvalue equation satisfied by the spherical har- 
monics. The vector and tensor harmonics of odd par- 



l rn 
A 



-efD B Y lm and X 



rim 

AB 



Wa°D 1 



OAB vim 
^AB 



0. 



ity are X 

e§D A )D c Y lm - X% also is tracefree 

We first express the electric-type tidal potentials in 
terms of the even-parity spherical harmonics. We begin 
with £"\ which we decompose as 



£W(v,6 A ) =J2£$(v)Y lm (6 A ), 



(2.8) 



in terms of harmonic components £m (v). There are 2Z+1 
terms in the sum, and the 21 + 1 independent components 
of £l are in a one-to-one correspondence with the 21 + 1 

(0 



coeffici ents £ m ' ■ Returning to the original representation 
of Eq. (2.4 1, we find after differentiation that D A £^ — 
lil A £ a j~iil ~~ , 



A') 



and we conclude that 
1 



D A £ 



; / j °rn 1 j 



Ini 
A ■ 



(2.9) 



Em B ™ y ' m (^)' Differentiating the first expression, 
multiplying this by the L evi-C ivita tensor, and in- 



volving the second of Eqs. (|2.3| returns e A D B B^ = 



L-l 



From this we conclude that 



= y(-e/i? B )B« = ^B$>X!T. (2.11) 



A second differentiation yields —e A D B DcB"> = 
IcabB^ + 1(1 - l)n a A e apq nPB q bL „ 2 n B n L - 2 , and after 
symmetrization we obtain 



B 



(0 

AB 



1(1 

2 



1) 



(efD B +e B c D A )D c B^ 



AB- 



(2.12) 



III. EXTERNAL PROBLEM 

A. Even-parity sector 

In this subsection we determine the tidal deformation 
of the metric outside the matter distribution, in the even- 
parity sector. The unperturbed external solution is the 
Schwarzschild metric 



ds n 



-fdv 2 + 2dvdr + r 2 dVL 2 



(3.1) 



with / := 1 — 2M jr and M denoting the body's mass; 
the metric is valid for r > R, where R is the body's 
radius. We employ the perturbation formalism of Martel 
and Poisson [T7], and implement the light-cone gauge of 
Preston and Poisson [5]. 

In the light-cone gauge the even-parity metric pertur- 
bation is given by 



h l ™{r)Y lm {9% 



PvA 



PAB 



\ ~» A 



j l v m (r)Y A m (9 A ), 



^ 2 E 



K lm (r)n AB Y lm (9 A ) 



„2 \ " nl 



G Lm (r)Y A ^(6 A ). 



(3.2a) 
(3.2b) 

(3.2c) 



An addit ional differentiation using the last of 
Eqs. (O) reveals that D A D B £ m = -1£Iab£ {1) + 
1(1 — l)Q A fl B £ a b l-2^ L ~ 2 ■ From this we conclude that 

2 —[D A D B + h(l + l)n AB 



c(0 _ 

■'AB — 



1(1 



£ {l) 



1(1-1) 



- 1 AB- 



(2.10) 



We next express the magnetic-type potentials in 
terms of the odd-parity spherical harmonics. We be- 
gin with f?W := BiS^ L and its decomposition B"> — 



We consider each Z-mode separately, and we henceforth 
omit the label lui on the perturbation variables h vv , j v , 
K, and G, which depend on r only. As discussed by Pre- 
ston and Poisson, K lm can always be set equal to zero 
when the perturbation satisfies the vacuum field equa- 
tions; this represents a refinement of the light-cone gauge, 
and we shall make this choice here. 

To simplify the task of solving the field equations we 
set 



(1-1)1 



r l ei(r)££>, 



(3.3a) 



9 



Jv 

G 



(1-1)1(1 + 1) 
4 

(I- l)l 2 (l + l) 



e 4 (r)£ 



(0 

m ' 



r'e 7 (r)£. 



(0 

m i 



(3.3b) 
(3.3c) 



where the functions ei(r), e 4 (r), and e 7 (r) are to be de- 
termined. Substitution of Eqs. (3.3) into Eq. (3.2) pro- 
duces 



(1-1)1 
2 



l ex(r)£ 



(i) 



(3.4a) 



PvA = 



(l-l)(l+l) 

2 



r l+1 e 4 (r)£ A l} , (3.4b) 



PAB 



i(i + r 



r l+2 e 7 (r)E% 



(3.4c) 



where £ A , a nd are the tidal potentials intro- 
duced in Eq. (j^ . 

The motivation behind the introduction of the func- 
tions ei, e 4 , and e 7 goes as follows. We first observe that 
when we set e% = e 4 = e 7 = 1, the perturbation de- 



fined by Eqs. (3.3 1 or Eqs. (3.4 1 satisfies the equations of 



linearized theory for a perturbation of Minkowski space- 
time. This exercise reveals that h vv must be proportional 
to t 1 , j v to r 



i+i 



and G to r l ; the relative numerical coef- 



ficients between these fields are also determined by solv- 
ing the perturbation equations in flat spacetime. The 
remaining absolute numerical coefficient that relates the 
perturbation to the tidal moment £l is determined by 
the definition of the tidal moment in terms of the Weyl 
tensor of the perturbed spacetime; this coefficient — the 
factor — 2/[(Z — 1)1] in h vv — can be read off Eq. (3.26a) 
of Ref. [9]. 

Inserting the functions e\, e 4 , and e 7 in Eqs. (3.3) al- 



lows the perturbation to be a solution to the Einstein 
field equations linearized about the Schwarzschild metric 
instead of the Minkowski metric. We impose the bound- 
ary conditions 



ei(r — > oo) = e 4 (r 



30 



) = e? (r -*■ oo) = 1. (3.5) 



The field equations do not determine these functions 
uniquely. The light-cone gauge comes with a class of 
residual gauge transformations that preserve the light- 
cone nature of the coordinate system (see Preston and 
Poisson [8 ) . In the even-parity sector, and for static per- 
turbations, the residual gauge freedom that keeps K = 
is a one-parameter family described by 

ei -> ei - la 1 (2M/r) l+2 , (3.6a) 
e 4 e 4 + oi [(J - l)(l + 2) + 4M/r] (2M/r) l+1 , (3.6b) 
e 7 -> e 7 + 2Zai(2M/r) m , (3.6c) 

in which a\ is the (dimensionless) parameter. The resid- 
ual gauge freedom does not interfere with the boundary 
conditions of Eq. (3.51. 



When K is allowed to change, the residual gauge free- 
dom becomes a three-parameter family. In this case we 



have 

ei -> e x - Zai(2M/r) /+2 + a 3 (2M/r) t+2 , (3.7a) 
e 4 -> e 4 + ai [(Z - 1)(Z + 2) + 4M/r] (2M/r) ;+1 

-(i + l)a 3 (2M/r) ,+1 , (3.7b) 
e 7 -> e 7 + 2Zai(2M/r) i+1 + 2a 2 (2M/r) 1 , (3.7c) 

and A' becomes 



K 



^[a 2 + a 3 (2M/r)]^. (3.8) 



Here a2 and 03 are two additional gauge parameters. 

The differential equations satisfied by e\, e 4 , and e 7 
can be extracted from the perturbation equations. These 
equations are coupled, and some effort must be devoted 
to their decoupling before an attempt is made to find 
solutions. We shall not describe these routine steps here. 
We state simply that the solutions are the ones that were 
displayed in Eqs. (1.6) and (1.7). These are given in a 



minimal implementation of the light-cone gauge, in which 
all constants of integrations are set equal to zero. The 
most general form of the solution is obtained from this by 
effecting the shifts described by Eqs. (3.7 1 and (3.8 1. The 
functions A n and B n are displayed for selected values of 
I in Table I. 

The metric perturbation can be represented in terms of 
gauge-invariant variables. We employ the set defined by 
Eqs. (4.10)-(4.12) of Martel and Poisson [If. According 
to these equations, and as can be directly verified from 
Eq. (13. 71), the variables 



2M 



h vv \ — h vv -\- ^ jv Af fG , 

h vr : MG' 

h rr := 2rG' + r 2 G", 



(3.9a) 

(3.9b) 
(3.9c) 



K := --j v + h(l + l)G + rfG' (3.96) 
r 2 

are gauge- invariant; a prime indicates differentiation with 
respect to r. We express them as 



K 



(I 


-1)1 


2 






1)1 




4 


(I 


-1)1 




2 


(I 


-1)1 



r e vv( r )£m 1 



r e rr (r)£^J , 



r l e K (r)£i\ 



(3.10a) 
(3.10b) 
(3.10c) 
(3.10d) 



in terms of new radial functions e vv , e vr , e rri and ex- 
Calculation reveals that these are given in terms of the 
old ones by 



1 2M 1 2M „ 

e 4 - i—r — M 



1 + 1 



/ + 1 r 



10 



TABLE I: Functions A n and B n for selected values of I, expressed in terms of z := 2M/r. The numbers fii and A; are given by 
A, = (21)1(21 + !)!/[(/ - 2)1(1 - !)!(/ + !)!(/ + 2)!] and m = (/ + l)Xi/l. 



1 = 

Ai 
A 4 
A 7 


2 


(I-,) 2 

1-2 

1-iz 2 


M2 = 30, A 2 = 20 

z 5 Bi = -jU 2 Ai ln(l - z) - f z(2 - z)(6 - 6z - z 2 ) 
z 5 B 4 = A 2 A 4 ln(l - z) + |z(12 - 6z - 2z 2 - z 3 ) 
z 5 B 7 = -M2A 7 ln(l - z) - 5z(6 + 3z - z 2 ) 


I = 
Ai 
A 4 
A 7 


3 


i(l-z) 2 (2-z) 
±(l-z)(3-2z) 
1-z+^z 3 


= 840, A 3 = 630 
z 7 B 1 - -ll 3 Ax ln(l - z) - 7z(120 - 240z + 130z 2 - 10z 3 - z 4 ) 
z 7 B 4 - A3A4 ln(l - z) + |z(180 - 210z + 30z 2 + 5z 3 + z 4 ) 
z 7 B 7 = -M3A 7 ln(l - z) - 14z(60 - 30z - 10z 2 + z 3 ) 


Ai 
A 4 
A 7 


4 


- z) 2 (14- 14z + 3z 2 ) 
i(l - a) (28 - 35z + 10z 2 ) 

± 3 z -(- 7 z 42 z 


^i 4 — 17 u40, A 4 — 14 112 

z 9 Bi = -At 4 Ai ln(l - z) - 21z(2 - z)(420 - 840z + 440z 2 - 20z 3 - z 4 ) 
z a B 4 - A4A4 ln(l - z) + 4 ^z(1680 - 2940z + 1370z 2 - 90z 3 - 9z 4 - z 5 ) 
z 9 B 7 - -m 4 A 7 ln(l - z) - 14z(1260 - 1470z + 270z 2 + 65z 3 - 3z 4 ) 


i = 
Ai 
A 4 
A 7 


5 


i(l - z) 2 (2 - z)(6 - 6z + z 2 ) 
55,(1 - z)(30 - 54z + 30z 2 - 5z 3 ) 
1 — £ z + 5 2 2 5_z 3 -l — z 5 

- 1 - 4 ^ T 3 z 12 ^ ^ 168 * 


Ms = 332 640, A 5 = 277 200 

z 11 B 1 = -//, 5 Ai ln(l - z) - 66z(5040 - 15120z + 16380z 2 - 7560z 3 + 1288z 4 - 28z 5 - z 6 ) 
z 11 B 4 = A5A4 ln(l - z) + 22z(12600 - 28980z + 21840z 2 - 5670z 3 + 210z 4 + 14z 5 + z 6 ) 
z al -B 7 = -fi 5 A 7 ln(l - z) - 66z(5040 - 8820z + 4410z 2 - 420z 3 - 77z 4 + 2z 5 ) 



1 



1(1 + 1) 
1 



2M/e' 7 , 

1 2M 



e 4 



e7 - 



i + 1 
1 



l + l r 



-e 7 



Z(Z- 



-2Me 7 , 



1 



2 // 

: r e 7> 



i + 1 



e 4 



i(i + l) 
-L(j + 3-4M/r)e 7 



Z(i+1) 



(3.11a) 

(3.11b) 
(3.11c) 

(3.11d) 



It is easy to see that these functions, like the old ones, 
all go to one as r goes to infinity. 

Substitution of our expressions for e±, 04, and e 7 into 
Eqs. (3.11 1 and repeated use of the properties of hyper- 



geometric functions reveal that 

e vv = fe vr = fe rr =A X + 2k cl (R/r) 2l+1 B 1 (3.12) 
and 

e K = A 7 + 2k cl (R/r) 2l+1 B 7 . (3.13) 

Notice that e vv , fe vr , and f 2 e rr are all equal to the mini- 
mal implementation of e l5 and Bk is equal to the minimal 
implementation of e 7 . All this shows that the relativistic 
Love numbers k c \ possess gauge-invariant significance. 



B. Odd-parity sector 

In the light-cone gauge the odd-parity metric pertur- 
bation is given by 



PvA = ^2h l ™(r)X l j n (9 A ), 



(3.14a) 



Pab 



(3.14b) 



We consider each Z-mode separately, and we henceforth 
omit the label Im on the perturbation variables h v and 
hi, which depend on r only. To simplify the task of 
solving the field equations we set 



3(1-1)1 
4 

3(1 -l)l 2 '' 



r l+ Mr)Bg\ 



J+2 



Hr)Bi 



(3.15a) 
(3.15b) 



where the functions 6 4 (r) and 67 (r) are to be determined. 
Substitution of Eqs. (3.151 into Eq. (3.141 produces 



2 

sj^iy 



PAB 



J+2 



where B^ and Bj^ B are the tidal potentials first intro- 



3/ 
3O 



b 7 (r)B 



(i) 

AB' 



(3.16a) 
(3.16b) 



duced in Eq. (2.7 1 



The motivation behind the introduction of the func- 
tions 64 and 67 is identical to what was done in the 
even-parity sector. When we set 64 = 67 = 1, the per- 
turbation defined by Eqs. (3.15) or Eqs. (3.16) satisfies 



the equations of linearized theory for a perturbation of 
Minkowski spacetime. This exercise reveals the relative 
numerical coefficients between h v and Y12 ■ The remaining 
absolute numerical coefficient that relates the perturba- 
tion to the tidal moment Bl is determined by the defi- 
nition of the tidal moment in terms of the Weyl tensor 
of the perturbed spacetime; this coefficient — the fac- 
tor 2/ [3(1 - 1)1] in h v — can be read off Eq. (3.26b) of 
Ref. 0. 



Inserting the functions 64 and 67 in Eqs. (3.151 allows 



the perturbation to be a solution to the Einstein field 
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equations linearized about the Schwarzschild metric in- 
stead of the Minkowski metric. We impose the boundary 
conditions 



&4(Y — * oo) = bi(r — ► oo) = 1. 



(3.17) 



The field equations do not determine these functions 
uniquely. As in the even-parity case we have a residual 
gauge freedom that preserves the nature of the light-cone 
coordinates. It is described by 



64 — > 64 
67 — * bj + a 



/2M\ l 



(3.18a) 
(3.18b) 



in which a is a (dimensionless) parameter. The resid- 
ual gauge freedom does not interfere with the boundary 
conditions of Eq. ( 3.17| ). 

The differential equations satisfied by 64 and 67 can be 
extracted from the perturbation equations. The solutions 
are displayed in Eqs. (1.6 1 and (1.7 1. They are given in a 



minimal implementation of the light-cone gauge, in which 
all constants of integrations are set equal to zero. The 
most general form of the solution is obtained from this 
by effecting the shifts described by Eqs. (3.18). 



The metric perturbation can be represented in terms 
of gauge-invariant variables. We employ the set defined 
by Eq. (5.7) of Martel and Poisson [17]. According to 
this, and as can be directly verified from Eq. (3.181, the 
variables 



h v . — h Vj 

h r := —hi — —/in 
r 2 

are gauge-invariant. We express them as 
2 



hr 



3(1 - 1)1 
2 



3(1-1)1 



r l+ %(r)B%, 
r i+1 6 r (r)£« 



(3.19a) 
(3.19b) 

(3.20a) 

(3.20b) 
(3.20c) 



in terms of new radial functions b v and b r . Calculation 
reveals that these are given in terms of the old ones by 



b v = 64, 
b r = b-; 



b'r- 



(3.21a) 
(3.21b) 



It is easy to see that these functions, like the old ones, 
all go to one as r goes to infinity. 

Substitution of our expressions for 64 and 67 into 



Eqs. (3.21 1 and repeated use of the properties of the hy- 



pergeometric functions reveal that 



by — fb r — A4 2 fcmaf 



(R/r) 2l+1 B 4 



(3.22) 



Notice that b v and fb r are both equal to 64, which is 
gauge-invariant. This shows that the relativistic Love 
numbers fc mag possess gauge-invariant significance. 



IV. INTERNAL PROBLEM 

A. Background metric for relativistic stellar 
models 

We begin with an examination of the internal gravi- 
tational field of a body that is not yet perturbed by an 
external tidal field. The body is spherically-symmetric, 
and the matter consists of a perfect fluid. In light-cone 
coordinates (v,r, 9 A ) the metric is expressed as 

dsl = -e 2l/ 7 dv 2 + 2e^ dvdr + r 2 dtt 2 , (4.1) 

with / = 1 — 2m(r)/r and ip = ip( r )- The Einstein field 
equations are 

m' = Anr 2 p : ip' = ~j~ (p + p) , (4.2) 

and the equation of hydrostatic equilibrium is 
m + 4irr 3 p 



P 



r 2 f 



-(P + P)- 



(4.3) 



Here p is the fluid's proper energy density, and p is the 
pressure. 

These equations can be integrated once an equation of 
state is specified. The boundary conditions are m(r = 
0) = and ip(r = 0) = i/jq, where ipo is chosen so that ip 
vanishes at the stellar surface: ip(r = R) = 0. 



B. Light-cone gauge 

The internal light-cone gauge is a modified version of 
the external gauge constructed by Preston and Poisson 
[8]. We define it properly in this section. 

The metric of Eq. (4.1 1 reveals the meaning of the coor- 
dinates (v,r,6 1 ) in the background spacetime. We note 
first that l a = —d a v is a null vector, so that the surfaces 
v — constant are null hypersurfaces; they describe light 
cones that converge toward r — 0. The vector 



r = (0,-e-*,0,0) 



(4.4) 



is tangent to the null generators of these light cones, and 
the expression reveals that 9 A is constant along the gen- 
erators. In addition, the affine parameter A that runs 
along the generators is related to r by dX = — dr. In 
the interior portion of the spacetime, r is no longer an 
affine parameter on the null generators; but it still pos- 
sesses the property of being an areal radius, in the sense 
that the area of a surface of constant (v,r) is given by 

47T7- 2 . 

In the internal light-cone gauge, the metric of the per- 
turbed spacetime is presented in coordinates (v,r,9 A ) 
that possess the same geometrical meaning as in the 
background spacetime. In particular, v continues to la- 
bel null hypersurfaces. 9 A continues to be constant along 
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the null generators, and r continues to be related to the 
affine parameter by dX = — dr. It is easy to show that 
these statements imply the same conditions 



Pvr 



Pr 



PrA = 



(4.5) 



that were employed in the external problem. The nonva- 
nishing components of the metric perturbation are there- 
fore p vv , p v A, and pab ■ The radial coordinate, however, 
will lose its meaning as an areal radius in the stellar in- 
terior. 

In the even-p arity sector the perturbation is decom- 
posed as in Eq. Q, and the fields h l ™, j l v m , K lm , G lm 
depend (in general) on the coordinates (v,r). An even- 
parity gauge transformation is generated by the vector 
field 3 a) with components 



rim ( (\A\ 



s„ = £e>,r)r< 

Ira 

* r = Y;£ n (v,r)Y lm {6 A ), 

Ira 

HA = £^,r)F|™(^). 



(4.6a) 
(4.6b) 
(4.6c) 



It can be shown that the condition h vr = determines 
£.„, that h rr = determines and that j r — deter- 
mines £. The gauge, however, is not determined uniquely 
There exists a residual gauge freedom that preserves the 
geometrical meaning of the coordinates. In the case of 
w-independent perturbations, the residual gauge freedom 
is a three-parameter family described by 



& = -a ie 2v 7 + a 2 , 
£ r = aie^, 

...2 



-air 



(4.7a) 
(4.7b) 

(4.7c) 



and £, as well 



Here we suppressed the Im labels on , £ 
as the constants a±, a,2, and a.3. 

In the odd-p arity sector the perturbation is decom- 
posed as in Eq. (3.141, and the fields h l ™, h ™ depend (in 
general) on the coordinates (v,r). An odd-parity gauge 
transformation is generated by the vector field E a , with 
components 



Z r = 0, 



(4.8) 



It can be shown that the condition h r = determines 
£. In this case also there exists a residual gauge free- 
dom that preserves the geometrical meaning of the co- 
ordinates. In the case of ^-independent perturbations, 
the residual gauge freedom is a one-parameter family de- 
scribed by 



ar 



(4.9) 



Here also we suppressed the Im labels on £ and the con- 
stant a. 



3.2 1 and (3.14) can be used 



The decompositions of Eq 
to compute 8G a f3, the perturbation of the Einstein ten- 
sor inside the body. The even-parity sector decouples 
from the odd-parity sector, and the perturbation takes 
the form of 



Irn 
Irn 

EA lm\rlm 
rr 1 5 

Im 
Irn 



SG rr 
5G v a 

SG r A 



Im 
A 



Im 
A 



(4.10a) 
(4.10b) 
(4.10c) 
(4.10d) 
(4.10e) 



°v A A J> 

SG AB = J2( A i m ^A B Y lm + A\ m Y AB + B lm X%). 

(4.10f) 



im 



Here the even-parity fields A vv , A vr , A rr , A V1 A ri A\,, 
A$ and the odd-parity fields B v , B r , B depend on v and 
r only. In the case of a stationary perturbation, they 
depend on r only. 

The expressions are too long to be displayed here. 
In practice they are easily generated with GRTensorll 
|18j by specializing the perturbation to an axisym- 
metric mode m = with a specific multipole order 
I, With Y lm = Y{6) we have Yg = Y' , Y^ = 0, 
Y ge = -cos0y'/sin0 - + 1)Y, Y e<p = 0, and 
= sin 6 cos &Y' + |Z(i + 1) sin 2 9Y in the even-parity 
case, and Xg = 0, X^ = sin 9Y' , Xgg = 0, X 0( p — 
-cos9Y' - + 1) sin AY", and X^ = in the odd- 
parity case. The definition of the metric implements 
the constraint Y" = ~ cos 6Y' / s'm 6 - 1(1 + l)Y on the 
spherical-harmonic functions, and this simplifies the final 
expression for the perturbed Einstein tensor. 



C. Energy- momentum tensor 

We consider stationary tides raised by a tidal environ- 
ment characterized by an electric-type tidal moment El 
and a magnetic- type tidal moment Bl', these are actu- 
ally time-dependent, but the dependence is sufficiently 
slow that it can be neglected in the process of integrat- 
ing the Einstein field equations. The perturbed metric 
will therefore carry a parametric dependence upon v. 

The fluid's velocity vector in the background configu- 
ration is given by u a = (e - ^/ -1 / 2 , 0, 0, 0). In the per- 
turbed configuration it becomes u a = (u v , 0, 0, 0), re- 
flecting the fact that the tide is stationary and does not 
create motion within the fluid. The time component of 
the vector changes by virtue of the fact that the met- 
ric changes; we have that u v 
Su v = \t 



i e -3V f-3/2 



-4>f-i/2 + Su v^ with 
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After lowering the index on u a with the perturbed 
metric g^p + p a p, we find that u v = —e^f 1 / 2 + 8u v , 
u r = /~ 1//2 + 8u r , and ua = Sua, with 



Su v = 


1 -V f-V2„ 


(4.11a) 


Su r = 


i -2f f -3/2 

2^ J i^DHi 


(4.11b) 


Su A = 


e^T^PvA- 


(4.11c) 



These expressions are valid in the light-cone gauge. The 
perturbation Sua can be decomposed into even-parity 
and odd-parity components; the perturbations Su v and 
Su r are necessarily of even parity. 

The perturbation in the energy-momentum tensor is 
generated by the perturbation in u a , but also by a per- 
turbation in the density p and pressure p created by the 
tide; these are related by the equation of state. We have 

ST a/3 = (p+p)(u a Su /3 + U l3 Su a ) +pp a /3 

+ (Sp + 5p)u a Uf3 + {Sp)g al3 , (4.12) 



and in the light-cone gauge this reads 



ST VV 



-PPv 



2V 



ftp, 



ST vr — — Sp, 
ST vA = -PPvA, 

ST„ = (p + p)e-^T 2 Pw + r\Sp + Sp), 
ST rA = e-^f- 1 (p + p)p vA , 
STab = PPab + r 2 6ptt A B- 



(4.13a) 
(4.13b) 
(4.13c) 
(4.13d) 
(4.13e) 
(4.13f) 



The perturbations ST v a, ST t a, and STab can be decom- 
posed into even-parity and odd-parity components; the 
perturbations 5T VV , ST vr , and ST rr are necessarily of even 
parity. 

From Eqs. ( |4.13 1 we find that ST a p is given by 



ST VV 
ST VT 

ST vA 
5T rA 
ST A B 



■I III 



rim 



I III 

J2qI?y 

I III 

J2q i ™y 

1 rn 
in 1 

£(<#"%' 

lm 



A 



rlrr 
■ A 



plm j^lm\ 



(4.14a) 
(4.14b) 
(4.14c) 
(4.14d) 
(4.14c) 



q -y 



AB 



AB 



in 1 



The even-parity fields are 



Qw = -ph vv + e 2v> / a 



-e* a, 



(4.14f) 

(4.15a) 
(4.15b) 





= (p + p)e- 2 ^'f- 2 h vv + f-\cr + q), 


(4.15c) 




= -pjv, 


(4.15d) 


Q r 


= e-*>r 1 (p + p) Jv , 


(4.15e) 




= r 2 (pK + q), 


(4.15f) 




= r 2 pG, 


(4.15g) 



and the perturbations in the density and pressure were 
also decomposed in spherical harmonics: 



Sp = J2 a lm Y lm , Sp=J2 q lm Y lm 

lm lm 

The odd-parity fields are 

Py - - pfl V ) 

P r = e-^r 1 (p + p) h v , 
P = ph 2 . 



(4.16) 



(4.17a) 
(4.17b) 
(4.17c) 



Information about Sp and Sp, or a and q, can be 
obtained from the equation of hydrostatic equilibrium. 
In the perturbed spacetime the equation states that 
(p + p)a a + d a p — 0, where p = p + Sp is the perturbed 
density p = p + Sp is the perturbed pressure, and a a 
is the perturbed acceleration of the fluid elements. The 
equation becomes 



(p + p)5a a + (Sp + Sp)a a + d a Sp = 



(4.18) 



when expressed in terms of the perturbations Sp, Sp, 
and Sa a . The unperturbed acceleration has a r = 
i e -2-4ij-i^ e -4ify ^ j£ g on }y nonvanishing component, 

and the perturbation has components 

5a v = 0, (4.19a) 



Sa r 



Sqa 



1 



-2ip 



f 1 d A p v 



(4.19b) 
(4.19c) 



Substitution of Eqs. ( 4.16[) a nd (4.191, as well as p v 
E/™ h l ™Y lm , into Eq. (|4.18| reveals that 



and 



-\(p + p)e-^r 2 (e 2 ^fYh l 
-le-^f-'^fYio + q) 



-(p + p)e- 2 ^f- 1 h vv . 



(4.20) 
(4.21) 



If we next differentiate Eq. ( |4.21 1 and insert the result 
within Eq. ( |4.20 1, we discover that 



(p + pyh vv = -(e^fY(<j + q) 



(4.22) 



The last two equations allow us to express a lm and q lm 
directly in terms of h 1 ™; hydrostatic equilibrium implies 
that these are not independent variables. 
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D. Perturbation equations: even-parity sector 

The useful combinations of Einstein field equations are 

Ex := (A vv - 8ttQ vv ) + e^f(A vr - 8nQ vr ) = 0, 

(4.23a) 

E 2 := {A vv - 8nQ vv ) + 2e lp f(A vr - 8irQ vr ) 

+ e^f (A rr - 8nQ rr ) = 0, (4.23b) 

E 3 :=(A rr - 8nQ rr ) = 0, (4.23c) 

E 4 := e-^rE 2 + 2f(A v - 8ttQ v ) 

+ 2e i 'f (A r - 8nQ r ) = 0. (4.23d) 

These are a set of coupled differential equations for the 
variables h vv (r), j v (r), K(r), and G(r); the remaining 
field equations are redundant by virtue of the Bianchi 
identities. The explicit forms reveal that Ex = is a first- 
order differential equation for j v , E 2 — is a first-order 
differential equation for h vv , E 3 — is a second-order 
differential equation for K, and E4 = is a first-order 
differential equation for G. 

The field equations can be manipulated to yield a de- 
coupled equation for the master function 

h vv := h vv + e-*(e^f)'j v -\r*f{eWf)'G' 



2e^ 



(m + 47rr 3 p)j tI 
e 2,( 'f(m + ATrr 3 p)G'. 



(4.24) 



This function is gauge invariant, and it joins smoothly 
with the external version of Eq. (3.9 1 at r = R. The 
master equation is 



where 



A = 

B = 



r 2 K v + Arh' vv - Bh vv = 0, (4.25) 



1 _ _ 27rr 2 (p + 3p) , (4.26a) 
r 

l(l + l)-A-Kr 2 {p + p)(?,+ d p) . (4.26b) 



The master equation is equivalent to Eq. (27) of Ref. [T^] . 
in which H := e~ 2 ^ f~ 1 h vv is used as an alternative 
choice of dependent variable. 

The master equation can be derived by the following 
procedure. First, integrate the field equation E$ := A$ — 
8irQ$ = and obtain j v = ^r 2 fe^G'. This implies that 
h vv = h vv . Second, make the substitution in the other 
field equations. The result is that Ex now involves h vv , 
G', and G"\ E 2 involves h vv , h' vv , K, K' ', and G'\ E 3 
involves h vv , K', and K"; and Ex involves h vv , K, K' , 
G, and G' . Third, differentiate E 2 with respect to r, and 
use Ex to eliminate the terms in G" , and E 3 to eliminate 
the terms in K" . The result is that E' 2 now involves h vv , 



h' vv , h" v , K, K' , G, and G' . Fourth, construct the linear 
combination rE' 2 +aE 2 +bEx and determine the functions 
a and b that eliminate all terms involving K, K', G, G' . 
The solution is unique, and the final result is Eq. (4.251. 



For numerical integration it is advantageous to make 
the same substitution as in Eq. (3.101, 



(1-1)1 



re. 



(0 



(4.27) 



and to rewrite Eq. ( |4.25[ ) as a second-order differential 
equation for e vv (r). This function joins smoothly with 
the external version of Eq. (3.12), and k c \ is determined 



by matching the values of the internal and external func- 
tions (along with their first derivatives) at r = R. 



E. Perturbation equations: odd-parity sector 



The useful combinations of field equations are 



Ox 

o 2 



(B v - 8irP v ) 
(B v - 8irP v ) 



0, 

e^f(B r 



87rP r ) 



(4.28a) 
0. (4.28b) 



The first is a second-order differential equation for h v , 
while the second is a first-order differential equation for 
h 2 . 

The equation Ox = is fully decoupled, and the per- 
turbation variable h v is easily shown to be gauge invari- 
ant, as it was in the external problem. The master vari- 
able for the odd-parity sector is therefore h v :— h v , and 
the master equation is 



" 2 h" 



Frh' - Gh v = 0, 



where 



F 



G = 



zl 7 7 7 

1(1 + 1) + 8irr 2 (p + p) 

r 



(4.29) 

(4.30a) 
(4.30b) 



This equation is equivalent to Eq. (31) of Ref. [12], in 
which ip := rh! v — 2h v is used as an alternative choice of 
dependent variable. The function h v joins smoothly with 
the external version of Eq. (3.19 1 at r = R. 



For numerical integration it is advantageous to make 
the same substitution as in Eq. (3.201, 



J+ 



3(1-1)1 



(4.31) 



and to rewrite Eq. ( |4.29[ ) as a second-order differential 
equation for b v (r). This function joins smoothly with the 
external version of Eq. (3.221, and fc mag is determined by 
matching the values of the internal and external functions 
(along with their first derivatives) at r = R. 
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IMPLEMENTATION FOR POLYTROPES 



The relativistic Love numbers k e \ and A; mag are deter- 



mined by the numerical integration of Eqs. (4.251 and 
(4.29 1 and matching with the external solutions at r = R. 



This defines a simple computational procedure that can 
be implemented for any choice of equation of state. In 
this section we describe the steps that are involved when 
the polytropic form 



p = Kp 



l + l/n 



(5.1) 



is adopted; here K and the polytropic index n are con- 
stants. We choose, however, to deviate from the pro- 
cedure just outlined: Instead of integrating the master 
equations for the variables h vv and h v , we integrate the 
complete set of independent field equations. This allows 
us to calculate all components of the metric perturbation, 
and matching them across r = R determines, in addition 
to the Love numbers, the gauge parameters oi, «2 5 ^3, 
and a that are automatically selected by the internal so- 
lution. 1 



p c = b n+1 /K n , and 6 can be used in place of p c to label a 
stellar model, given a choice (K, n) of equation of state. 
We also note the relation m /r = (n+l)b. We find that 
the units mo and tq vary with b even when the equation 
of state is fixed. To eliminate this dependence it is useful 
to define the alternative units 



M 



("+l) 3 V /2 

\/47r 



Ro 



n + 1 
4n 



K n ' 2 , (5.5) 



which do not depend on b. We have that mo — 

Mob (3-n)/2 and r() = # fc(l-n)/2. 

In terms of the dimensionless variables the field equa- 



tions (4.2 1 and (4.3 1 become 



2 an 



dp 

— = (n + l)6 , 

d£ f 

d9 _ _{p + b^9 n+1 ){l + b9) 



(5.6a) 
(5.6b) 
(5.6c) 



A. Unperturbed stellar model 



The numerical integration of Eqs. (4.2 1 and (4.3 1 is 



conveniently accomplished by introducing the dimension- 
less variables 9, p, and £ defined by 



P = Pc 



P = p c n+ \ 



= mop,, r = r £. (5.2) 



Here p c :— p(r — 0) is the central density, and p c :— 
Kpc +1 ^ n is the central pressure. The units of mass and 
radius are chosen to be 



m := 47rr p c , 



(n + l)p c 
47rp2 



(5.3) 



so as to simplify the form of the field equations. 
It is useful to introduce also a "relativistic factor" 



b :=p c /pc, 



(5.4) 



which determines the degree with which the stellar model 
is relativistic. In terms of this we have p c = b n /K n , 

1 There is no strong rationale for proceeding in this way. The 
honest truth is that we became aware of Eq. | |4.25| only after 
completing the numerical work. We derived the master equation 



with / = 1 — 2(n+ l)bp/£. The boundary conditions are 
0(£ = 0) = 1, = 0) = 0, and = 0) = t^o- In 
the limit 6—^0 the model becomes nonrelativistic, and 
the equations for p and 9 can be combined into the well- 
known Lane-Emden equation; in the limit the equation 
for ip becomes irrelevant. 

The formulation of Eq. (5.6) is not optimal from a 
numerical point of view. For accurate integrations it is 
better to use the variable v := p/£, 3 instead of p, and 
x := ln£ instead of £. The system of equations becomes 



dv 

dx 
dtp 

dx 

de 

dx 



= 9 n ~ 3u, (5.7a) 
= (n+l)b£ 2 f- 1 6 n (l + b6), (5.7b) 
= -ff- X {v+b6 n+l ){l + b6), (5.7c) 



with / = 1 — 2(n + \)b^ 2 v. The integration begins at a 
large and negative value of x, so that £ = e x is small, 
with the starting values 

after noticing its appearance in Refs. [2] 112] and wondering why 
our formulation was more complicated than theirs. 



v = 3 - ^(1 + 6)(1 + 36)£ 2 + —(1 + 6)(1 + 36) [8n - 5 + (18n - 20)6 + (15 + 30n)6 2 ]£ 4 + 0(£ 6 ), (5.8a) 

9 = \-\{l + b)(l + 36)£ 2 + -^-(1 + 6)(1 + 36) [3n - 2nb + (30 + 15n)6 2 l£ 4 + 0(£ 6 ), (5.8b) 
6 360 

ifi = VA) + \(n + 1)6(1 + 6)£ 2 - i(n + 1)6(1 + 6) [n - 36 + (3 + 3n)6 2 ]£ 4 + O(f). (5.8c) 
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The integration stops at £ = £1, where 9 goes to zero, 
and ipo is chosen so that ip(^i) = 0. The stellar mass and 
radius are then given by 

M = M & (8-n)/a £?i/(£i), R = R b {1 - n)/2 £ u (5.9) 



in the units of Eq. ( |5.5| . The compactness of the body 
is measured by C := 2M/R = 2(n + l)6£ 2 ^(£i); this is 
dimcnsionless, and therefore independent of the units M 
and Rq. 



1)6/(4tit§), p c = (n+ l)& 2 /(47rr 2 ), and m = (n + l)br , 
with and v (as well as ip) depending on £. Finally, we 
use the fact that a term p' in the perturbation equations 
can be related to p' by the equation p' = {dp / dp)p' , with 
dp/ dp determined by the equation of state. 

Another useful set of substitutions is the one displayed 
in Eqs. (3.3 I, along with 



r$K = 



(Z-l)Z(Z + 2)(Z + 3) 



r l+2 e w (0£«\ (5.10) 



B. Perturbation: Even-parity sector 



The perturbation equations (4.23 1 are simplified by in- 



volving the background field equations (4.2 1 and (4.3 1 



They are also simplified by making the substitutions of 
Eqs. (5.2 1, (5.3 1, and (5.4 1; we therefore write p = p c 9 n , 
p = p c 8 n+1 , r = ro^, and m — toq£ 3 ^j where p c = (n + 



in which we replace the original variables with the radial 
functions e\, en, e 7 , and eio- These replacements are 
motivated by an analysis of the perturbation equations 
for small values of r, which reveals that h vv behaves as 
r l , j v as r +x , G as r , and K as r l+2 . The numerical 
factors in front of eio is inserted to simplify the form of 
the small-r expansion of K, as we shall see below. 
The final expression of the perturbation equations is 



= E 1 
= E 2 

= E 3 
= E 4 



-^ 4 + (l+l)e-^r 1 e 1 -r 1 A 1 e*, 

<e'i + lf~ l A 2 ex + le ,p f- 1 B 2 e 4 - 1(1 - l)(l + 2)e 2 ^ e 7 



-e 2 ^ 2 £ 2 e 10 



^B 2 f e 



3 „/ 



= -£e' 7 



2(Z + 3)~ " n ~ ±u (Z + 2)(Z + 3)~ ~" ~ 10 ' 
ee'{ + (Z + 2)(Z + 3)e- 2 ^r 2 A 3 e x - (I + 2)/- 1 B 3 e 10 - f^Cat e' 10 , 
Z(Z + 1) „_o . z 



2(Z- l)(Z + 2) 



e-^f- 2 A i e 1 + 



-l[l + 3-4(n+l)b^y]f- 1 e r + 



(Z-l)(Z + 2) 
1(1 + 1) 



e~^r 2 B^ 



1(1+1) 



2(1- l)(Z + 2)(Z + 3) 
[(n + l)b(v + bO n+1 )]f- i e eio, 



(I- l)(Z + 2) 2 (Z + 3) 

where a prime indicates differentiation with respect to £, and 

Z + 1 - 2(n + 1)Z>£ 2 [(Z + 2)v + b8 n+1 ], 
(I - 2)(l + 1) + 2(n + 1)6£ 2 [2(1 + l)u - Q n (l + 
1 - (n+ l)b^ 2 (u -b9 n+1 ), 
l-3 + 2(n + l)b£ 2 (v-b9 n+1 ), 



Ax 
A 2 
B 2 
C 2 
A 3 
B 3 
C 3 
Ai 
B 4 
C 4 



2/m+l 



nO 11 - 1 + (An + 3)b9 n + 3(n + 1)6 
Z + 3 - (n + 1)6£ 2 [2(1 + 3)i/ + 6 n (l + b9)} . 
2(1 + 3) - (n + 1)&£ 2 [4(Z + 3)i/ + 0" (1 + i 
(Z - 1)(Z + 2) + 2(n + 1)6£ 2 [2i/ - 9 n (l + , 
(Z - 1)(Z + 2) - (n + 1)Z>£ 2 [(Z 2 + Z - 4)i/ - Z(Z + l)66» n+1 ] 
Z-l-2(n+l)6£ 2 (z/ + &(9" +1 ). 



(5.11a) 

(5.11b) 
(5.11c) 



(5.11d) 



(5.12a) 
(5.12b) 
(5.12c) 
(5.12d) 
(5.12e) 
(5.12f) 
(5-12g) 
(5.12h) 
(5.12i) 
(5-12j) 



A small-£ expansion of these equations, using Eqs. (5. 
reveals that e\ = a + 0(£ 2 ), e 4 = a n e~^° 



o(e), e 7 



a e- 2 ^ +O(e), and e w = a e- 2 ^(l + 6)[3(n+l)6+n] + 
0(f 2 ), where ao is a parameter that must be determined 
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by matching the internal and external perturbations at 
the stellar boundary. 

The perturbation equations are easily written as a 
first-order dynamical system for the variables u\ := e±, 
112 '■= e4, u 3 = ey, U4 := eio, and U5 := £e' 10 . The nu- 
merical integration is carried out with x := ln£ as the 
independent variable, and the differential equations are 
integrated simultaneously with Eqs. (5.7 1 to determine 



the unperturbed stellar model. The integration proceeds 
from a large and negative value of x, for which £ = e x is 
small, and it stops at at £ = £1 where 9 goes to zero. 

The term n0 n_1 in A 3 originates from a term involving 
dp /dp oc 9~ l that multiplies p oc 9 n in the field equation 
for K (or eio). This term diverges at the stellar bound- 
ary when n < 1. The singularity is integrable, however, 
and it can be shown that the solution for K(r) (or eio) is 
actually well-behaved at the boundary. The divergence 
of ^3 nevertheless causes issues in the numerical integra- 
tion of the perturbation equations. For this reason, the 
accuracy achieved for n < 1 is limited compared with the 
accuracy obtained for n > 1. 

The internal perturbation must match the external 
perturbation at £ = £1 , or r — R, the position of the stel- 
lar boundary. The five internal functions ei, e 4 , e 7 , eio, 
and £e' 10 depend on one free parameter do- The external 
functions, on the other hand, depend on three gauge pa- 
rameters ai, a 2 , and 03, as well as the electric- type Love 
number k c \. The five matching conditions determine the 
five parameters uniquely, including the Love number. 

We suppose that the internal functions u 4l ■ ■ ■ ,u 5 are 
determined by setting ao = 1 in the numerical integra- 
tions. The desired functions ei, • • • , eio then differ from 
these by an overall multiplicative factor that we denote 
A -1 . We have 



e l 


= A~Vi, 


(5.13a) 


< 4 


= A _1 it 2 , 


(5.13b) 


e 7 


= A _1 it 3 , 


(5.13c) 


e 10 


= A _1 1t4, 


(5.13d) 


de% 
<% 


= A _1 w 5 , 


(5.13e) 



and the matching conditions are 



e 7 
e 10 

d£ 



out 

C 10 I 



e 10 



(5.14a) 
(5.14b) 
(5.14c) 
(5.14d) 

(5.14e) 



where each side of the equation is evaluated at £ = £1. 
The external expressions for ei, and e-j are presented 
in Eqs. (1.6 1 and (1.7 1, and these must be modified by 



the gauge adjustments of Eqs. (3.7 1 



The function e 10 is related to K by Eq. (5.101, and 



-10 



equation and its derivative with respect to r imply that 
at £ = £1, 

2(1 + 2)(Z + 3) C'^f 2 [02 + a 3 (2M/R)] , 

(5.15a) 

-2(/ + 2)(/ + 3)C'£f 2 [(; + 2)a 2 

+ (; + 3)a 3 (2M/i?)], (5.15b) 

where C := 2M/R is the compactness factor. These 
equations can be solved for a 2 and 03. Involving also the 
matching equations and Eqs. (5.131, we arrive at 



P 2 

\ ?i 

2 ~ 2(l + 2)(l + 3)C l 



(I + 3)u 4 + u 5 



(5.16a) 



Aa 3 = - 



2{l + 2)(l + 'S)C l + 1 



(l + 2) Ui + u 5 . (5.16b) 



We see that the gauge parameters a 2 and 03, rescaled by 
the unknown coefficient A, are determined by the numer- 
ical values obtained for it 4 and u 5 . 

To solve the remaining matching equations we transfer 
the a 2 and 03 terms from the right-hand side of Eqs. (3.7 1 



to the left-hand side. Taking Eqs. (5.161 into account, we 
form the combinations 



Ul 



w 2 := u 2 



w 3 



U3 



eg 


2(1 + 


2)(l + 


3) 


(I 


+m 




2(1 + 


m+ 


3) 




$ 


( 



(I + 2)u 4 + u 5 



(I + 2)u 4 + u 5 



(l + 2)(l + 3) 



(5.17a) 
(5.17b) 
(5.17c) 



which can be determined numerically. Involving now 
Eqs. ( 1.6 1, the matching conditions take the explicit form 



Wi = A 1 -X + 2Bf (Afc cl ) - IC ■ (XC l+1 ai ), (5.18a) 
w 2 = Ai ■ X - 2^—^Bi ■ (Afe e i) 

+ [(I - l)(l + 2) + 2C] • (AC i+1 ai), (5.18b) 
w 3 = A 7 ■ A + 2B 7 ■ (Afc c i) + 21 ■ (XC l+1 ai )- (5.18c) 

in these expressions the functions A n and B n are evalu- 
ated at r = R, or 2M/r = C. 

If we define a vector w — (wi,W2,w 3 ) of numerical 
quantities, and another vector p = (A, Afc re i, AC' +1 ai) of 
unknown parameters, these equations take the form of 
the matrix equation w = M p, with a matrix M that is 
known analytically. Solving for p, the Love number is 
finally determined by k e \ = p 2 /pi . 



C. Perturbation: Odd-parity sector 

To arrive at the final form of the perturbation equa- 
tions (4.281 we follow the same steps as in the even- 



the external expression for K is given by Eq. (3.8 1. This parity sector. These include making the substitutions 
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of Eqs. (3.151, to replace the original variables h v and h 2 
with the radial functions 64 and 67. 
The perturbation equations are 

= Oi = -£ 2 6 4 ' - b' A + f^Gi 64, (5.19a) 

= 2 = -£6' 7 -lb 7 + le-^T 1 64, (5.19b) 

with 

F l = 2(1 + 1) - (n + 1)6£ 2 [4(Z + l)y + n (l + b9)} , 

(5.20a) 

d = (n + l)bf [2(1 - l)(l + 2)v + (l + 3)9 n (l + bQ)} . 

(5.20b) 

A small-£ expansion of these equations reveals that 64 — 
a +O(£ 2 ) and 67 — aoe~^°+0(£ 2 ), where ao is a param- 
eter that must be determined by matching the internal 
and external perturbations at the stellar boundary. 

The perturbation equations are easily written as a first- 
order dynamical system for the variables v± := 64, v 2 := 

£6' 4 , and v 3 '■= b 7- 

The internal perturbation must match the external 
perturbation at £ = £1, or r = R, the position of the 
stellar boundary. The three internal functions 64, £b' A , 
and 67 depend on one free parameter ao. The exter- 
nal functions, on the other hand, depend on one gauge 
parameter a as well as the magnetic-type Love number 
fc mag . The three matching conditions determine the three 
parameters uniquely, including the Love number. 

We suppose that the perturbation equations for v\, v 2 , 
and v$ are integrated with ao = 1. The desired internal 
functions 64 and 67 are then given by 



bf = A~V, 
V? = A -1 u 3 , 



(5.21a) 
(5.21b) 
(5.21c) 



where A is an unknown constant. The matching condi- 
tions are 



tout 
°4 : 



.dbf 



db° 4 ut 



b™ = b° u 



(5.22a) 
(5.22b) 
(5.22c) 



where each side of the equation is evaluated at £ = £1. 
The external expressions for 64 and 67 are presented 



in Eqs. (1.6), together with the gauge adjustment of 
Eq. (3.18). We observe that b± is gauge-invariant, and 



that the purpose of the matching equation for b-j is to 
determine the (uninteresting) gauge parameter a. 

We focus on the two equations involving 6 4 . Using 
Eqs. ( 1.6 1, we find that the explicit form of the matching 



conditions is 



Vi 



Aa-X 



2 — - — Ba ■ (Afc mag ), 



v 2 = 



CA' A ■ A + 2^— [CB' A + (21 + l)B 4 ] ■ (Xk n 



(5.23b) 



In these expressions the functions Aa, A' a := dA^/dz, B4, 
and B' 4 := dB^jdz are evaluated at z :— 2M/r = C. 

If we define a vector v — (vi,v 2 ) of numerical quan- 
tities, and another vector p — (A, A/c mag ) of unknown 
parameters, these equations take the form of the ma- 
trix equation v = Mp, with a matrix M that is known 
analytically. Solving for p, the Love number is finally 
determined by fc mag = p 2 /pi- 

To evaluate the derivatives of A 4 and B4 with respect 
to z we use the well-known property of hypergeometric 
functions that (d/dz)F(a,b; c; z) — (ab/c)F(a + 1,6 + 
l;c+ l;z). 



VI. NUMERICAL RESULTS 

The computations presented in this work were gener- 
ated with two independent codes, one written by each au- 
thor. Consistency between our results provides evidence 
that each set of computations were carried our correctly, 
and the comparison allows us to estimate the numerical 
accuracy of our results. 

The background spacetime is constructed by solving 
the Einstein field equations for a spherical matter config- 
uration with a polytropic equation of state. The equa- 
tions were formulated in Sec. V A, and the system of 



equations (5.7 1 is integrated numerically for selected val- 
ues of the polytropic index n. The integration begins at 
a large and negative value of the radial variable x = ln£, 
using the starting values listed in Eqs. (5.8 1. It proceeds 



until 6 changes sign at the stellar boundary, x = X\. 
In the first code, the integration is performed using the 
Bulirsh-Stoer method as implemented in the Numeri- 
cal Recipes routine bsstep, which is embedded within 
odeint; we use the Second Edition of Numerical Recipes 
[T9"] , and the code is written in C++. In the second code, 
the integration is performed using the embedded Rungc- 
Kutta Prince-Dormand method as implemented in the 
GNU Scientific Library routine rk8pd, which is embed- 
ded within odeiv; we use version 1.9 of the libraries [2"U] . 
and the code is written in C. In each code all floating- 
point operations are carried out with double precision. 
The accuracy of the integration is determined by the in- 
tegrator's tolerance e and the errors of order £ 6 that are 
incorporated in the starting values. As Eqs. (5.7) are 



(5.23a) 



exceptionally well conditioned toward numerical integra- 
tion, a high degree of accuracy can easily be achieved. 
We estimate that our stellar configurations are computed 
accurately to at least twelve significant digits. 

The stellar boundary is identified with the help of a 
bisection search for the solution to 9(x) = 0. In the first 
code this is carried out with the Numerical Recipe rou- 
tine zbrent; the search is loosely bracketed between the 
values x < x\ (where 9 is positive) and x 2 > x\ (where 
9 is negative). In the second code this is carried out with 
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the GNU Scientific Library routine brent, using a similar 
bracketing method. The search is carried out with high 
accuracy, again of the order of twelve significant digits. 



The even-parity perturbation equations (5.11 ) are next 
integrated for selected values of n and I, simultaneously 
with the background field equations (5.7 1. Once more 
the integration begins at a large and negative value of 
x, using the starting values derived in Sec. V B, and it 
proceeds up to x = x \. In the first code we continue to 
use bsstep and odeint, and the caption of Table II dis- 
cusses the accuracy of these integrations. In the second 
code we continue to use rk8pd and odeiv; the tolerance 
of the integrator is set uniformly to e = 1.0e-12, and all 
integrations begin at x = —10.0. Each code returns the 
values of u\, U2, U3, 114, and U5 at the stellar boundary. 



The odd-parity equations (5.191 are integrated in ex- 
actly the same way. Here the codes return the values of 
Vi, t>2, and V3 at the stellar boundary. 



The matching problem of Eqs. (5.181 requires the nu- 

l\p, where 
is known 



merical solution of the matrix equation w 
w is constructed from the perturbations, 
analytically, and p is the vector of unknown parame- 
ters, which include the electric-type Love number k e \. 
In the first code the system of equations is solved by 
performing an LU decomposition of the matrix M, and 
this is handled by the Numerical Recipes routines ludcmp 
and lubksb. In the second code the LU decomposi- 
tion is handled by the GNU Scientific Library routines 
gsl_linalg_LU_decomp and gsl_linalg_LU_solve. In 
view of the small number of equations involved (three) , 
this task is essentially carried out at machine precision. 
The final output is k c \. 



TABLE II: Integration errors for even-parity perturbations. 
For each selected value of n the first row shows the value of e, 
the integrator's tolerance. When e = 1.0e-12 the integrations 
are started at x = —7.0, so that the errors in the starting 
values are of the order of 1.0e-12. When e > 1.0e-12 the in- 
tegrations are started at x = —6.5, so that the errors in the 
starting values are of the order of 1.0e-ll. For the odd-parity 
equations the tolerance of the integrator is set uniformly to 
e = 1.0e-12, and all integrations begin at x = —7.0. The 
second column shows 8, an intrinsic measure of the accuracy 
of our results. This is defined as S := |f mo dci — ^pcrtl/fmodei, 
where iv dci i s the value of v at the stellar boundary £ = £1 as 
determined with exquisite precision by integrating the stellar- 
model equations only, while f por t is the value as determined 
by also integrating the perturbation equations. The least ac- 
curate determinations are for small values of b; the accuracy 
typically improves by two orders of magnitude at larger values 
of b. For reasons that were explained at the end of Sec. IV D, 
when n < 1 the accuracy that can be achieved for the even- 
parity perturbations is more limited than what is achieved for 
v\ for these cases S gives an overestimate of the true accuracy. 
For n > 1, and for the odd-parity perturbations, S should be 
an accurate measure of our accuracy. 
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The matching problem of Eqs. (5.231 is handled in ex- 



actly the same way. Here the final output is the magnetic- 
type Love number fc mag . 



Our results are presented in the figures displayed in 
Sec. I and in the tables provided in the Appendix. The 
electric-type and magnetic-type Love numbers are com- 
puted for selected values of n and I, as functions of the 
relativistic parameter b := p c / p c and the compactness 
C := 2M/R. The allowed interval begins at b = and 
C = 0, where the equations reduce to their Newtonian 
limit, and ends at b — 6 max and C — C max , where the 
stellar configuration achieves its maximum mass. Each 
table caption discusses the estimated accuracy of our re- 
sults. Overall we claim an approximate accuracy of nine 
significant digits for the Love numbers (with some excep- 
tions, as detailed in the table captions). 
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APPENDIX: TABLES OF RELATIVISTIC LOVE 
NUMBERS 



TABLE III: Love numbers for n = 0.50 and I = 2. Inte- 
gration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fed = 4.491539995415e-01. This provides evidence that 
our results for the electric-type Love numbers are accurate 
to five significant digits. We believe that our results for the 
magnetic-type Love numbers are accurate to nine significant 
digits. 



b 


2M/R 


feel 






0.0000000000 


0.0000000000 


4.4915295584c- 


■01 


O.OOOOOOOOOOc+00 


0.0162962963 


0.0627859865 


3.6857599573e- 


•01 


1.6896831556c-03 


0.0651851852 


0.2085406132 


2.2117103643e- 


•01 


4.1596525372c-03 


0.1466666667 


0.3636165454 


1.1528493484e- 


•01 


4.8522931593e-03 


0.2607407407 


0.4883414066 


6.0195686393e- 


-02 


4.2583493395c-03 


0.4074074074 


0.5772867923 


3.3825660571e- 


-02 


3.3759103748c-03 


0.5866666667 


0.6379537736 


2.0879830949e- 


•02 


2.6271666539e-03 


0.7985185185 


0.6789539591 


1.4107892752e- 


-02 


2.0809271632c-03 


1.0429629630 


0.7068171264 


1.0311605798e- 


-02 


1.7004161103e-03 


1.3200000000 


0.7259502382 


8.0453742292c- 


■03 


1.4372173622e-03 



TABLE V: Love numbers for n = 1.00 and I = 2. Inte- 
gration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fc c i = 2.599088771480e-01. This provides evidence that 
our results for the electric-type Love numbers are accurate 
to nine significant digits. We believe that our results for the 
magnetic-type Love numbers are also accurate to nine signif- 
icant digits. 
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0.0000000000 


0.0000000000 
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O.OOOOOOOOOOc+00 


0.0054320988 
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■01 
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•01 
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•01 


2.3418562643e-03 


0.0869135802 


0.2449940757 


1.0135201470e- 


•01 


2.9092025152c-03 


0.1358024691 


0.3264977638 


6.8656738911e- 


•02 


3.0470778168c-03 


0.1955555556 


0.3971100356 


4.6672564713c- 
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0.4400000000 
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1. 7360105151c- 


■02 


1.9854445481c-03 



TABLE IV: Love numbers for n = 0.75 and I = 2. Inte- 
gration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fed = 3.434291771770e-01. This provides evidence that 
our results for the electric-type Love numbers are accurate 
to nine significant digits. We believe that our results for the 
magnetic-type Love numbers are also accurate to nine signif- 
icant digits. 



0.0000000000 
0.0092592593 
0.0370370370 
0.0833333333 
0.1481481481 
0.2314814815 
0.3333333333 
0.4537037037 
0.5925925926 
0.7500000000 



2M/R 

0.0000000000 
0.0363293144 
0.1294393332 
0.2455878442 
0.3564603549 
0.4485200887 
0.5194393410 
0.5719839827 
0.6101589410 
0.6376107260 



feci 

3.4342917761e-01 
3.0528456672c-01 
2.2114677984e-01 
1.4055702918c-01 
8.4890119064c-02 
5.1708383719c-02 
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TABLE VI: Love numbers for n = 1.25 and / = 2. Inte- 
gration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns k e \ = 1.943393766752e-01. This provides evidence that 
our results for the electric-type Love numbers are accurate 
to nine significant digits. We believe that our results for the 
magnetic-type Love numbers are also accurate to nine signif- 
icant digits. 
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-02 


1.8682931012c-03 
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TABLE VII: Love numbers for n = 1.50 and I = 2. Inte- 
gration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fc c i = 1.432787706403e-01. This provides evidence that 
our results for the electric-type Love numbers are accurate 
to nine significant digits. We believe that our results for the 
magnetic-type Love numbers are also accurate to nine signif- 
icant digits. 
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TABLE IX: Love numbers for n = 2.00 and / = 2. Inte- 
gration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fc c i = 7.393839192094e-02. This provides evidence that 
our results for the electric-type Love numbers are accurate 
to nine significant digits. We believe that our results for the 
magnetic-type Love numbers are also accurate to nine signif- 
icant digits. 
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TABLE VIII: Love numbers for n = 1.75 and I = 2. Inte- 
gration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fe c i = 1.039154459896e-01. This provides evidence that 
our results for the electric-type Love numbers are accurate 
to nine significant digits. We believe that our results for the 
magnetic-type Love numbers are also accurate to nine signif- 
icant digits. 



0.0000000000 
0.0018518519 
0.0074074074 
0.0166666667 
0.0296296296 
0.0462962963 
0.0666666667 
0.0907407407 
0.1185185185 
0.1500000000 



2M/R 

0.0000000000 
0.0064743298 
0.0251788386 
0.0541259986 
0.0904962838 
0.1311832131 
0.1732771212 
0.2143825635 
0.2527489948 
0.2872548584 



feci 

1.0391544596e-01 
1.0123582505c-01 
9.3755924816e-02 
8.2923108986e-02 
7.0531274857c-02 
5.8177787311c-02 
4.6952247447e-02 
3.7393783155e-02 
2.9616069949e-02 
2.3478510254e-02 



"■mag 

O.OOOOOOOOOOe+00 

6.7882541528c-05 

2.5279292539e-04 

5.0667193661c-04 

7.7147592064c-04 

9.9862752282c-04 

1.1599057338c-03 

1.2480829694e-03 

1.2710091094e-03 

1.2440214762e-03 



TABLE X: Love numbers for n = 0.50 and I = 3. Inte- 
gration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fed = 2.033844048605e-01. This provides evidence that 
our results for the electric-type Love numbers are accurate 
to five significant digits. We believe that our results for the 
magnetic-type Love numbers are accurate to nine significant 
digits. 



b 


2M/R 


feci 




femag 


0.0000000000 


0.0000000000 


2.0338399420e- 


■01 


0.0000000000e+00 


0.0162962963 


0.0627859865 


1.5613095764e- 


■01 


7.6806695868c-04 


0.0651851852 


0.2085406132 


7.9298498872e- 


■02 


1.5334802180e-03 


0.1466666667 


0.3636165454 


3.3876635518e- 


-02 


1.4053484595c-03 


0.2607407407 


0.4883414066 


1.4803140981& 


-02 


1.0040207606c-03 


0.4074074074 


0.5772867923 


7.2690559784c- 


-03 


6.8628867282c-04 


0.5866666667 


0.6379537736 


4.0970043659e- 


■03 


4.8495681685c-04 


0.7985185185 


0.6789539591 


2.6182028430e- 


■03 


3.6226354828e-04 


1.0429629630 


0.7068171264 


1.8555610194& 


■03 


2.8623307831c-04 


1.3200000000 


0.7259502382 


1. 4266248200c- 


■03 


2.3757226341c-04 
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TABLE XI: Love numbers for n = 0.75 and I = 3. Integration 
of the Newtonian Clairaut equation (for 6 = 0) returns fc e i = 
1.479565910794e-01, and this value was copied in the first row 
of the Table. [We were not able to accurately compute the 
electric-type Love number for 6 = for these specific values 
of n and I. The reason has to do with the fact that for these 
values, £e' 10 = 0(£ 4 ) instead of being of order £ 2 near £ = 0; 
the integrator then has difficulty moving out of the small- 
£ region and the number of steps required exceeds the set 
limit.] We believe that our results for the electric-type Love 
numbers are accurate to nine significant digits, and that our 
results for the magnetic-type Love numbers are also accurate 
to nine significant digits. 



2M/R 



0.0000000000 
0.0092592593 
0.0370370370 
0.0833333333 
0.1481481481 
0.2314814815 
0.3333333333 
0.4537037037 
0.5925925926 
0.7500000000 



0.0000000000 
0.0363293144 
0.1294393332 
0.2455878442 
0.3564603549 
0.4485200887 
0.5194393410 
0.5719839827 
0.6101589410 
0.6376107260 



1.4795659108c-01 
1.2656912474c-01 
8.2729118342e-02 
4.5802954407e-02 
2.3980661301o-02 
1.2870072782e-02 
7.3958196272o-03 
4.6219154497e-03 
3. 1388009813o-03 
2.2971603620c-03 



O.OOOOOOOOOOc+00 
3.7606006609c-04 
9.7442657049e-04 
1. 1847104807o-03 
1. 0533665372o-03 
8. 2451036345o-04 
6.2355209034o-04 
4.7654009722o-04 
3. 7510397112o-04 
3.0589573877e-04 



TABLE XIII: Love numbers for n = 1.25 and I = 3. Inte- 
gration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns k cl = 7.558993098406e-02. This provides evidence that 
our results for the electric-type Love numbers are accurate to 
eight significant digits. We believe that our results for the 
magnetic-type Love numbers are accurate to nine significant 
digits. 



b 




2M/R 






k 


mag 





.0000000000 


0.0000000000 


7.5589930713c- 


02 





.OOOOOOOOOOc+00 





.0037037037 


0.0140910881 


7.0841000307e- 


-02 


9 


.1484438895c-05 





.0148148148 


0.0535218477 


5.8788456530e- 


■02 


3 


.02205752130-04 





.0333333333 


0.1109782733 


4. 4167684362c- 


-02 


5 


.05539401660-04 





.0592592593 


0.1774670027 


3.0965107128e- 


■02 


6 


.18486094290-04 





.0925925926 


0.2449789236 


2.0910331154c- 


-02 


6 


.34022707240-04 





.1333333333 


0.3079051472 


1.3985888466e- 


■02 


5 


.8644352242c-04 





.1814814815 


0.3631764781 


9.4660656773c- 


■03 


5. 


.12646004260-04 





.2370370370 


0.4096934360 


6.5780656202e- 


03 


4. 


3641911298e-04 





.3000000000 


0.4475972736 


4.7331767581c- 


03 


3 


.6881598982o-04 



TABLE XII: Love numbers for n = 1.00 and I = 3. Inte- 
gration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fed = 1.064540469774e-01. This provides evidence that 
our results for the electric-type Love numbers are accurate 
to nine significant digits. We believe that our results for the 
magnetic-type Love numbers are also accurate to nine signif- 
icant digits. 



0.0000000000 
0.0054320988 
0.0217283951 
0.0488888889 
0.0869135802 
0.1358024691 
0.1955555556 
0.2661728395 
0.3476543210 
0.4400000000 



2M/R 

0.0000000000 
0.0211887760 
0.0788326459 
0.1586178173 
0.2449940757 
0.3264977638 
0.3971100356 
0.4550360296 
0.5008905693 
0.5363092473 



feci 

1.0645404707e-01 
9.6920315090e-02 
7.4354157385e-02 
5.0164622016e-02 
3.1412814100e-02 
1. 9197799362o-02 
1.1894296947c-02 
7.6520308738e-03 
5.1742095921c-03 
3.6915952107e-03 



"■mag 

O.OOOOOOOOOOe+00 

1.7712436105c-04 

5.4046973625e-04 

8.0927197289e-04 

8.7805794150c-04 

8.0536598779c-04 

6.7961 132670c-04 

5.5466172884c-04 

4.5069370183e-04 

3.7043176559e-04 



TABLE XIV: Love numbers for n = 1.50 and / = 3. Inte- 
gration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fed = 5.284852444148e-02. This provides evidence that 
our results for the electric-type Love numbers are accurate to 
eight significant digits. We believe that our results for the 
magnetic-type Love numbers are accurate to nine significant 
digits. 



0.0000000000 
0.0025925926 
0.0103703704 
0.0233333333 
0.0414814815 
0.0648148148 
0.0933333333 
0.1270370370 
0.1659259259 
0.2100000000 



2M/R 

0.0000000000 
0.0095061309 
0.0366144717 
0.0775387088 
0.1272141307 
0.1805262693 
0.2332124569 
0.2822627839 
0.3259042474 
0.3633567807 



feci 

5.2848524127e-02 
5.0478328434o-02 
4.4172676367e-02 
3.5833784527e-02 
2.7409928071c-02 
2.0154523547e-02 
1.4516418814c-02 
1.0411153801e-02 
7.5326749926e-03 
5.5503018177e-03 



*-mag 

0.0000000000e+00 

4.6818602757c-05 

1.6314544473e-04 

2.9516806482e-04 

3.9553177536c-04 

4.4448482790c-04 

4.4702358096c-04 

4.1916508131c-04 

3.7680447176c-04 

3.3107508747e-04 
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TABLE XV: Love numbers for n = 1.75 and / = 3. Inte- 
gration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fc c i = 3.628620386492e-02. This provides evidence that 
our results for the electric-type Love numbers are accurate 
to nine significant digits. We believe that our results for the 
magnetic-type Love numbers are also accurate to nine signif- 
icant digits. 



2M/R 



0.0000000000 
0.0018518519 
0.0074074074 
0.0166666667 
0.0296296296 
0.0462962963 
0.0666666667 
0.0907407407 
0.1185185185 
0.1500000000 



0.0000000000 
0.0064743298 
0.0251788386 
0.0541259986 
0.0904962838 
0.1311832131 
0.1732771212 
0.2143825635 
0.2527489948 
0.2872548584 



3.6286203851c-02 
3.5106196263e-02 
3.1863068261e-02 
2. 7304450667o-02 
2.2304606782e-02 
1.7570120169e-02 
1. 3509063981o-02 
1.0255357455c-02 
7. 7654175911o-03 
5.9143273246c-03 



O.OOOOOOOOOOc+OO 

2.3637643620e-05 

8.5448967297e-05 

1.6345757744c-04 

2. 3443830054o-04 

2.83381 16595e-04 

3.0609684764o-04 

3.0629942762e-04 

2.9105710804c-04 

2.6737047585e-04 



TABLE XVII: Love numbers for n = 0.50 and I = 4. Inte- 
gration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns k cl = 1.250625809919e-01. This provides evidence that 
our results for the electric-type Love numbers are accurate 
to six significant digits. We believe that our results for the 
magnetic-type Love numbers are accurate to nine significant 
digits. 



2M/R 



0.0000000000 
0.0162962963 
0.0651851852 
0.1466666667 
0.2607407407 
0.4074074074 
0.5866666667 
0.7985185185 
1.0429629630 
1.3200000000 



0.0000000000 
0.0627859865 
0.2085406132 
0.3636165454 
0.4883414066 
0.5772867923 
0.6379537736 
0.6789539591 
0.7068171264 
0.7259502382 



1. 2506232752o-01 
8.9880099035e-02 
3.8670819375c-02 
1.3511973834e-02 
4.9038796705e-03 
2.0751263151o-03 
1.0476762526e-03 
6.1942259351c-04 
4.1610976670e-04 
3.0850496903c-04 



0.0000000000e+00 

4.1259713417c-04 

6. 8285490951o-04 

4.9948942099e-04 

2.9045885185e-04 

1.6878093386o-04 

1. 0600793048o-04 

7.2980383110O-05 

5.4565032992e-05 

4. 3643093032o-05 



TABLE XVI: Love numbers for n = 2.00 and I = 3. Inte- 
gration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fed = 2.439399851849e-02. This provides evidence that 
our results for the electric-type Love numbers are accurate 
to nine significant digits. We believe that our results for the 
magnetic-type Love numbers are also accurate to nine signif- 
icant digits. 



0.0000000000 
0.0013580247 
0.0054320988 
0.0122222222 
0.0217283951 
0.0339506173 
0.0488888889 
0.0665432099 
0.0869135802 
0.1100000000 



2M/R 

0.0000000000 
0.0044806121 
0.0175403971 
0.0381005822 
0.0645673693 
0.0950758472 
0.1277344111 
0.1608206440 
0.1929044035 
0.2228982181 



feci 

2.4393998521e-02 
2.3804409363o-02 
2.2147566853e-02 
1. 9717948676o-02 
1. 6892301431o-02 
1. 4022630411O-02 
1.1366324623e-02 
9.0665499267e-03 
7.1697064326e-03 
5.6581213050e-03 



"■mag 

O.OOOOOOOOOOe+00 

1.1834173013c-05 

4. 3858418141o-05 

8.7207170874c-05 

1.3133222716c-04 

1.6767660885c-04 

1.9163622562c-04 

2.0252062134c-04 

2.0229022040e-04 

1.9405485402e-04 



TABLE XVIII: Love numbers for n = 0.75 and I = 4. In- 
tegration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fed = 8.731859904775e-02. This provides evidence that 
our results for the electric-type Love numbers are accurate to 
eight significant digits. We believe that our results for the 
magnetic-type Love numbers are accurate to nine significant 
digits. 



0.0000000000 
0.0092592593 
0.0370370370 
0.0833333333 
0.1481481481 
0.2314814815 
0.3333333333 
0.4537037037 
0.5925925926 
0.7500000000 



2M/R 

0.0000000000 
0.0363293144 
0.1294393332 
0.2455878442 
0.3564603549 
0.4485200887 
0.5194393410 
0.5719839827 
0.6101589410 
0.6376107260 



feci 

8.7318599147e-02 
7.1916151340e-02 
4.2448736557e-02 
2.0467648174c-02 
9.2636229108e-03 
4.3517890992e-03 
2.2374880023c-03 
1.2812234382e-03 
8.1468581955c-04 
5.6833824831e-04 



^mag 

0.0000000000e+00 

1.9820540551o-04 

4.5682358629e-04 

4.7418211732e-04 

3. 5718798449o-04 

2.4045653422c-04 

1.6046670110O-04 

1.1123202446e-04 

8.1394894614c-05 

6.2954932476e-05 
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TABLE XIX: Love numbers for n = 1.00 and I = 4. Inte- 
gration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fc c i = 6.024125532418e-02. This provides evidence that 
our results for the electric-type Love numbers are accurate 
to nine significant digits. We believe that our results for the 
magnetic-type Love numbers are also accurate to nine signif- 
icant digits. 



6 2M/R kci fe. 



0.0000000000 





.0000000000 


6.0241255395c-02 


0. 


OOOOOOOOOOo+OO 


0.0054320988 





.0211887760 


5. 3646913671o-02 


9 


0169742307o-05 


0.0217283951 





.0788326459 


3.8686876721e-02 


2 


5599084572c-04 


0.0488888889 





.1586178173 


2.3850649236e-02 


3 


4511609518C-04 


0.0869135802 





.2449940757 


1.3456371026e-02 


3 


3171170926e-04 


0.1358024691 





.3264977638 


7.3994944976o-03 


2. 


.69138316580-04 


0.1955555556 





.3971100356 


4.1569144571o-03 


2 


.0272875503o-04 


0.2661728395 





.4550360296 


2.4561332329e-03 


1 


.4990535283c-04 


0.3476543210 





.5008905693 


1. 5478883581o-03 


1. 


,1224767948c-04 


0.4400000000 





.5363092473 


1.0442398816e-03 


8. 


.6440104267c-05 



TABLE XXI: Love numbers for n = 1.50 and I = 4. Inte- 
gration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fc e i = 2.739306738271e-02. This provides evidence that 
our results for the electric-type Love numbers are accurate to 
eight significant digits. We believe that our results for the 
magnetic-type Love numbers are accurate to nine significant 
digits. 



b 




2M/R 


k 


si 




&mag 





.0000000000 


0.0000000000 


2 


.7393067294c- 


02 


O.OOOOOOOOOOc+00 





.0025925926 


0.0095061309 


2 


.5904136864c- 


-02 


2.1590905746e-05 





.0103703704 


0.0366144717 


2 


.2022860127e- 


•02 


7.2661900733e-05 





.0233333333 


0.0775387088 


1 


,7085103298c- 


-02 


1.2456848415c-04 





.0414814815 


0.1272141307 


1 


2356917233e- 


•02 


1.5600796417c-04 





.0648148148 


0.1805262693 


8 


5360694947e- 


■03 


1.6255526538c-04 





.0933333333 


0.2332124569 


5 


7656653837e- 


•03 


1.5121350651e-04 





.1270370370 


0.2822627839 


3 


.8858798086e- 


■03 


1.3142526538c-04 





.1659259259 


0.3259042474 


2 


.6550241954e- 


03 


1.1010882007c-04 





.2100000000 


0.3633567807 


1 


.8598572889c- 


■03 


9.0858200272e-05 



TABLE XX: Love numbers for n = 1.25 and I = 4. Inte- 
gration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fed = 4.096746123839e-02. This provides evidence that 
our results for the electric-type Love numbers are accurate to 
eight significant digits. We believe that our results for the 
magnetic-type Love numbers are accurate to nine significant 
digits. 



b 


2M/R 


feel 




^mag 


0.0000000000 


0.0000000000 


4.0967461 120e- 


02 


O.OOOOOOOOOOe+00 


0.0037037037 


0.0140910881 


3.7831532658& 


•02 


4.4448609345e-05 


0.0148148148 


0.0535218477 


3.0102008525c- 


•02 


1.3969326868e-04 


0.0333333333 


0.1109782733 


2.1224824306e- 


-02 


2.1676154291c-04 


0.0592592593 


0.1774670027 


1.3778277197c- 


•02 


2.4208766237c-04 


0.0925925926 


0.2449789236 


8.5688261330e- 


•03 


2.2508402313c-04 


0.1333333333 


0.3079051472 


5.2857027758c- 


-03 


1.8910615297c-04 


0.1814814815 


0.3631764781 


3.3202258257e- 


•03 


1.5123549667c-04 


0.2370370370 


0.4096934360 


2.1611153793& 


-03 


1.1903913573c-04 


0.3000000000 


0.4475972736 


1.4718380690o-03 


9.4141547749e-05 



TABLE XXII: Love numbers for n = 1.75 and I = 4. Inte- 
gration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fed = 1.795919608352e-02. This provides evidence that 
our results for the electric-type Love numbers are accurate to 
eight significant digits. We believe that our results for the 
magnetic-type Love numbers are accurate to nine significant 
digits. 



b 


2M/R 






femag 


0.0000000000 


0.0000000000 


1.7959195798c- 


-02 


0.0000000000e+00 


0.0018518519 


0.0064743298 


1.7256414891e- 


-02 


1.0299492462c-05 


0.0074074074 


0.0251788386 


1. 5352372303c- 


•02 


3.6329038871c-05 


0.0166666667 


0.0541259986 


1.2748788134c- 


■02 


6.6858496307c-05 


0.0296296296 


0.0904962838 


1.0001909468e- 


■02 


9.1232683035c-05 


0.0462962963 


0.1311832131 


7.5212893015e- 


03 


1.0412456086c-04 


0.0666666667 


0.1732771212 


5.5036193453c- 


03 


1.0576629570c-04 


0.0907407407 


0.2143825635 


3.9751203447e- 


■03 


9.9452555069e-05 


0.1185185185 


0.2527489948 


2.8695308743e- 


03 


8.8982268992e-05 


0.1500000000 


0.2872548584 


2.0913685076e- 


■03 


7.7283983050e-05 
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TABLE XXIII: Love numbers for n = 2.00 and I = 4. In- 
tegration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fe c i = 1.150774963254e-02. This provides evidence that 
our results for the electric-type Love numbers are accurate 
to nine significant digits. We believe that our results for the 
magnetic-type Love numbers are also accurate to nine signif- 
icant digits. 



b 


2M/R 


feel 




&mag 


0.0000000000 


0.0000000000 


1.1507749634c- 


-02 


O.OOOOOOOOOOc+00 


0.0013580247 


0.0044806121 


1.1175986242c- 


-02 


4.8517346178e-06 


0.0054320988 


0.0175403971 


1.0253199166e- 


■02 


1.7666877492e-05 


0.0122222222 


0.0381005822 


8.9267030374e- 


■03 


3.4155212233e-05 


0.0217283951 


0.0645673693 


7.4272319821e- 


■03 


4.9576781437e-05 


0.0339506173 


0.0950758472 


5.9573506314c- 


■03 


6.0605384193c-05 


0.0488888889 


0.1277344111 


4.6508432777c- 


■03 


6.6035419880c-05 


0.0665432099 


0.1608206440 


3.5682657166c- 


■03 


6.6396107161c-05 


0.0869135802 


0.1929044035 


2.7150400670c- 


■03 


6.3100849384e-05 


0.1100000000 


0.2228982181 


2.0653404463c- 


■03 


5.7696936776e-05 



TABLE XXV: Love numbers for n = 0.75 and I = 5. Inte- 
gration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fc c i = 5.904211079675e-02. This provides evidence that 
our results for the electric-type Love numbers are accurate 
to nine significant digits. We believe that our results for the 
magnetic-type Love numbers are also accurate to nine signif- 
icant digits. 



b 


2M/R 


feci 




femag 


0.0000000000 


0.0000000000 


5.9042110830c- 


02 


O.OOOOOOOOOOc+00 


0.0092592593 


0.0363293144 


4.6830270753e- 


-02 


1.1653809273e-04 


0.0370370370 


0.1294393332 


2.4976158762e- 


•02 


2.4051097139e-04 


0.0833333333 


0.2455878442 


1. 0491664768c- 


•02 


2.1488760292e-04 


0.1481481481 


0.3564603549 


4.1030598266e- 


■03 


1.3812306450c-04 


0.2314814815 


0.4485200887 


1.6845383073e- 


•03 


8.0340315047c-05 


0.3333333333 


0.5194393410 


7.7279612957c- 


-04 


4.7382520827c-05 


0.4537037037 


0.5719839827 


4.0394347591c- 


■04 


2.9755158176c-05 


0.5925925926 


0.6101589410 


2.3943047563e- 


■04 


2.0179431020e-05 


0.7500000000 


0.6376107260 


1. 5846629836c- 


■04 


1.4743433223e-05 



TABLE XXIV: Love numbers for n = 0.50 and / = 5. In- 
tegration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fed = 8.758378097872e-02. This provides evidence that 
our results for the electric-type Love numbers are accurate 
to five significant digits. We believe that our results for the 
magnetic-type Love numbers are accurate to nine significant 
digits. 



0.0000000000 
0.0162962963 
0.0651851852 
0.1466666667 
0.2607407407 
0.4074074074 
0.5866666667 
0.7985185185 
1.0429629630 
1.3200000000 



2M/R 

0.0000000000 
0.0627859865 
0.2085406132 
0.3636165454 
0.4883414066 
0.5772867923 
0.6379537736 
0.6789539591 
0.7068171264 
0.7259502382 



feci 

8.7583597477e-02 
5.8953726923c-02 
2.1502135233e-02 
6.1438638016c-03 
1.8484475878c-03 
6.7146633825c-04 
3.0201770900c-04 
1.6415053975e-04 
1.0383693027c-04 
7.3777009618e-05 



O.OOOOOOOOOOc+00 

2.4566625412c-04 

3.4017649539e-04 

2.0040145035c-04 

9.5309205703c-05 

4.7055290188c-05 

2.6143552436c-05 

1.6470651693e-05 

1.1563913004e-05 

8.8490411630e-06 



TABLE XXVI: Love numbers for n = 1.00 and / = 5. In- 
tegration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fed = 3.929250022713e-02. This provides evidence that 
our results for the electric-type Love numbers are accurate 
to nine significant digits. We believe that our results for the 
magnetic-type Love numbers are also accurate to nine signif- 
icant digits. 



0.0000000000 
0.0054320988 
0.0217283951 
0.0488888889 
0.0869135802 
0.1358024691 
0.1955555556 
0.2661728395 
0.3476543210 
0.4400000000 



2M/R 

0.0000000000 
0.0211887760 
0.0788326459 
0.1586178173 
0.2449940757 
0.3264977638 
0.3971100356 
0.4550360296 
0.5008905693 
0.5363092473 



feci 

3.9292500283e-02 
3.4232186798e-02 
2.3214946910e-02 
1.3083821677c-02 
6.6517658359e-03 
3.2897043734c-03 
1.6738586645c-03 
9.0663862820e-04 
5.3118742797c-04 
3.3782253209e-04 



*-mag 

0.0000000000e+00 

5.1595279179c-05 

1.3691347868e-04 

1.6725789364c-04 

1.4339451212c-04 

1.0355791045c-04 

6.996284441 lc-05 

4.7012138969e-05 

3.2481406675c-05 

2.3432821346e-05 
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TABLE XXVII: Love numbers for n= 1.25 and / = 5. In- 
tegration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fc c i = 2.574776897544e-02. This provides evidence that 
our results for the electric-type Love numbers are accurate 
to nine significant digits. We believe that our results for the 
magnetic-type Love numbers are also accurate to nine signif- 
icant digits. 



b 


2M/R 


feel 




&mag 


0.0000000000 


0.0000000000 


2.5747768892c- 


■02 


O.OOOOOOOOOOc+00 


0.0037037037 


0.0140910881 


2.3431732289c 


■02 


2.4484594391e-05 


0.0148148148 


0.0535218477 


1.7882343045c 


■02 


7.3469913671c-05 


0.0333333333 


0.1109782733 


1.1838067300c- 


-02 


1.0630025814c-04 


0.0592592593 


0.1774670027 


7.1175929856c 


■03 


1.0903958368c-04 


0.0925925926 


0.2449789236 


4.0764785256c 


■03 


9.2520830998c-05 


0.1333333333 


0.3079051472 


2.3180122367c 


■03 


7.1006099307c-05 


0.1814814815 


0.3631764781 


1.3500205542c 


■03 


5.2196651860c-05 


0.2370370370 


0.4096934360 


8.2182992059c- 


-04 


3.8121417411e-05 


0.3000000000 


0.4475972736 


5.2873946327c- 


-04 


2.8280143688e-05 



TABLE XXX: Love numbers for n = 2.00 and I = 5. Inte- 
gration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns k e \ = 6.419966834096e-03. This provides evidence that 
our results for the electric-type Love numbers are accurate 
to nine significant digits. We believe that our results for the 
magnetic-type Love numbers are also accurate to nine signif- 
icant digits. 



b 


2M/R 


feci 




femag 


0.0000000000 


0.0000000000 


6.4199668350c- 


03 


O.OOOOOOOOOOc+00 


0.0013580247 


0.0044806121 


6.2054683813c 


■03 


2.3273253180e-06 


0.0054320988 


0.0175403971 


5. 6146736603c- 


•03 


8.3410581918e-06 


0.0122222222 


0.0381005822 


4.7814250834c 


-03 


1.5722362491e-05 


0.0217283951 


0.0645673693 


3.8647546556c 


-03 


2.2075823575c-05 


0.0339506173 


0.0950758472 


2.9960238568c 


■03 


2.5950234182c-05 


0.0488888889 


0.1277344111 


2.2531557240c 


■03 


2.7084179213c-05 


0.0665432099 


0.1608206440 


1.6628187985c 


03 


2.6037432070c-05 


0.0869135802 


0.1929044035 


1.2172373651c 


03 


2.3660856250e-05 


0.1100000000 


0.2228982181 


8.9227889200c- 


■04 


2.0720823774e-05 



TABLE XXVIII: Love numbers for n = 1.50 and I = 5. In- 
tegration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fe c i = 1.656876321404e-02. This provides evidence that 
our results for the electric-type Love numbers are accurate to 
eight significant digits. We believe that our results for the 
magnetic-type Love numbers are accurate to nine significant 
digits. 



b 


2M/R 


feci 




^mag 


0.0000000000 


0.0000000000 


1.6568763135c- 


-02 


O.OOOOOOOOOOe+00 


0.0025925926 


0.0095061309 


1.5513747923c 


-02 


1.1396769321e-05 


0.0103703704 


0.0366144717 


1.2817102827c 


-02 


3.7144252136e-05 


0.0233333333 


0.0775387088 


9.5123354123c- 


-03 


6.0593529815c-05 


0.0414814815 


0.1272141307 


6.5071687257c 


•03 


7.1292599566c-05 


0.0648148148 


0.1805262693 


4.2236890233c 


•03 


6.9268626197c-05 


0.0933333333 


0.2332124569 


2.6751245641c 


■03 


5.9938429831c-05 


0.1270370370 


0.2822627839 


1.6935347048c 


•03 


4.8543637402e-05 


0.1659259259 


0.3259042474 


1.0918134805c 


•03 


3.8082191743e-05 


0.2100000000 


0.3633567807 


7.2625290353c- 


-04 


2.9627611402e-05 



TABLE XXIX: Love numbers for n = 1.75 and I = 5. In- 
tegration of the Newtonian Clairaut equation (for 6 = 0) re- 
turns fe c i = 1.043995446810e-02. This provides evidence that 
our results for the electric-type Love numbers are accurate 
to nine significant digits. We believe that our results for the 
magnetic-type Love numbers are also accurate to nine signif- 
icant digits. 



b 


2M/R 


feci 




^mag 


0.0000000000 


0.0000000000 


1.0439954387c- 


-02 


O.OOOOOOOOOOc+00 


0.0018518519 


0.0064743298 


9.9635882364c 


•03 


5.1907364491c-06 


0.0074074074 


0.0251788386 


8.6906171854c- 


-03 


1.7903644602c-05 


0.0166666667 


0.0541259986 


6.9954658202c 


■03 


3.1806852682e-05 


0.0296296296 


0.0904962838 


5.2724767209c 


■03 


4.1473311159e-05 


0.0462962963 


0.1311832131 


3.7857836403c 


■03 


4.4914541070e-05 


0.0666666667 


0.1732771212 


2.6366800282c 


■03 


4.3129384841c-05 


0.0907407407 


0.2143825635 


1.8116921836c 


-03 


3.8310103373c-05 


0.1185185185 


0.2527489948 


1.2463184219c 


•03 


3.2437461311c-05 


0.1500000000 


0.2872548584 


8.6864853073c 


-04 


2.6760732768e-05 
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